Statistical mechanics of two-dimensional Euler flows and minimum enstrophy states 



A. Naso, P.H. Chavanis and B. Dubrulle 

^ Laboratoire de Physique, 
Ecole Normale Superieure de Lyon and CNRS (UMR 5672), 
46 allee d'ltalie, 69007 Lyon, France 
^ Laboratoire de Physique Theonque (IRSAMC), 
CNRS and UPS, Umversite de Toulouse, 
F-31062 Toulouse, France 
^ SPEC/IRAMIS/CEA Saclay, and CNRS (URA 2464), 
91191 Gif-sur-Yvette Cedex, France 

(Dated: To be included later) 

A simplified thermodynamic approach of the incompressible 2D Euler equation is considered 
based on the conservation of energy, circulation and microscopic enstrophy. Statistical equilibrium 
states are obtained by maximizing the Miller-Robert-Sommeria (MRS) entropy under these sole 
constraints. We assume that these constraints are selected by properties of forcing and dissipation. 
We find that the vorticity fluctuations are Gaussian while the mean fiow is characterized by a 
linear Zu — tp relationship. Furthermore, we prove that the maximization of entropy at fixed energy, 
circulation and microscopic enstrophy is equivalent to the minimization of macroscopic enstrophy 
at fixed energy and circulation. This provides a justification of the minimum enstrophy principle 
from statistical mechanics when only the microscopic enstrophy is conserved among the infinite class 
of Casimir constraints. Relaxation equations towards the statistical equilibrium state are derived. 
These equations can serve as numerical algorithms to determine maximum entropy or minimum 
enstrophy states. We use these relaxation equations to study geometry induced phase transitions 
in rectangular domains. In particular, we illustrate with the relaxation equations the transition 
between monopoles and dipoles predicted by Chavanis & Sommeria [J. Fluid. Mech. 314, 267 
(1996)]. We take into account stable as well as metastable states and show that metastable states 
are robust and have negative specific heats. This is the first evidence of negative specific heats in 
that context. We also argue that saddle points of entropy can be long-lived and play a role in the 
dynamics because the system may not spontaneously generate the perturbations that destabilize 
them. 
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I. INTRODUCTION 

Two-dimensional turbulence has the striking prop- 
erty of organizing spontaneously into large-scale coherent 
structures. These coherent structures correspond tojets 
and vortices in geophysical and astrophysical flows [J [^l- 
They can be reproduced in numerical simulations and 
laboratory experiments [1, Q in either forced or unforced 
situations. These coherent structures involve the pres- 
ence of a mean flow and fluctuations around it. The 
mean flow turns out to be a steady state of the pure 2D 
Euler equations (without forcing and dissipation) at some 
coarse-grained scale. In several cases, the steady state is 
characterized by a linear relationship between vorticity 
UJ (or potential vorticity q = u + h in the presence of 
a topography h) and stream function ijj. For example, 
the case of a linear q —Jp relationship was considered 
early by Fofonoff (1954) [6| 0jS 0, simple model of oceanic 
circulation. These "Fofonoff flows" were found to emerge 
naturally from random initial conditions in numerical ex- 
periments of forced and unforced 2D turbulence [7l-[TT|. 
However, these results are not expected to be general. 



There exists many other cases in 2D turbulence where 
the q — tp relationship is not linear. The understanding 
and prediction of these quasi stationary states (QSS), in 
forced and unforced situations, is still a challenging prob- 
lem. Different approaches have been proposed to describe 
these QSSs. 

In the case of forced flows (for example oceanic flows 
experiencing a forcing by the wind and a dissipation), Ni- 
iler (1966) ^ and Marshall & Nurser (1986) [3 have 
proposed that forcing and dissipation could equilibrate 
each other in average and determine a QSS that is a 
steady state of the unforced and inviscid 2D Euler equa- 
tion. This QSS is characterized by a functional relation- 
ship q = /('0) between potential vorticity and stream 
function where the function / is selected by the proper- 
ties of forcing and dissipation. In particular, these au- 
thors discussed the properties that forcing must possess 
to generate Fofonoff flows. 

In the case of unforced flows, a phenomenological ap- 
proach, called the minimum enstrophy pri ncip le, was pro- 
posed by Bretherton & Haidvogel (1976) [1J|. It is based 
on the inverse cascade [ssj process of Batchelor (1969) 
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[iSj. It is argued that, in the presence of a small vis- 
cosity, the (potential) enstrophy decays while the energy 
is approximately conserved. This is what Matthaeus & 
Montgomery (1980) [l^ have called selective decay. In 
that case, we can expect that the system will relax to- 
wards a state that minimizes (potential) enstrophy at 
fixed energy. This is mainly a postulate. This principle 
leads to a steady state of the 2D Euler equation charac- 
terized by a linear q — ip relationship. When applied to 
geophysical flows, this principle can give an alternative 
justification of Fofonoff flows. The minimum enstrophy 
principle has been generalized by Leith (1984) [l3| so as 
to take into account the conservation of angular momen- 
tum in order to describe isolated vortices. However, this 
principle is purely phenomenological and it is difficult to 
establish its domain of validity. 

A statistical mechanics approach of 2D turbulence has 
been developed by Kraichnan (1967,1975) [3 [3 based 
on the truncated Euler equations. In that case, the dy- 
namics conserves only the energy and the enstrophy (the 
other constraints of the Euler equation arc lost by the 
truncation). Since this system is Liouvillian, we can ap- 
ply the methods of statistical mechanics. This is the 
so-called energy-enstrophy statistical theory. In the ab- 
sence of topography, the truncated statistical mechanics 
predicts a homogeneous flow with an equilibrium energy 
spectrum of the form E{k) = k/(a + bk"^) correspond- 
ing to an equipartition distribution. In the presence of 
a topography (or /3-effect), this statistical mechanics ap- 
proach was generalized by Salmon, HoUoway & Hender- 
shott (1976) [l^]. It leads to a mean flow where the aver- 
age potential vorticity q and the average stream function 
Tp are related to each other by a linear relationship. When 
applied to geophysical flows, this approach provides a 
justification of Fofonoff flows from statistical mechanics. 
However, the fact that the truncated statistical mechan- 
ics breaks the conservation of some integrals of the 2D 
Euler equations may be considered as a limitation of this 
approach. 

Another statistical mechanics of 2D turbulence, based 
on a point vortex approximation, was initiated by On- 
sager (1949) |21| and further developed by Joyce & Mont- 
gomery (1973)[23l and Lundgren & Pointin (1977) [H] in 
a mean field approximation. The statistical equilibrium 
state maximizes the usual Boltzmann entropy (adapted 
to point vortices) while conserving the energy and the 
number of vortices of each species. This statistical me- 
chanics predicts an equilibrium state where the relation- 
ship between vorticity and stream function is given by 
a Boltzmann distribution, or a superposition of Boltz- 
mann distributions (on the different species). However, 
the point vortex approximation is a crude approximation 
of real turbulent flows where the vorticity field is contin- 
uous. 

A statistical theory of 2D turbulence valid for contin- 
uous vorticity fields has been develope d by Miller (1990) 
H] and Robert & Sommeria (1991) [25|. This theory 
takes into account all the constraints of the 2D Euler 



equation. The statistical equilibrium state maximizes a 
mixing entropy while conserving energy, circulation and 
all the Casimirs. This leads to equilibrium states with 
more general mean fiows than in the previous approaches. 
In particular, the oJ — -0 relationship and the fluctuations 
around it are determined by the initial conditions and 
can take various shapes. However, some connections with 
the earlier works can be found. For example, Robert & 
Sommeria (1991) [25| note that the results of the point 
vortex approach can be recovered in the dilute limit of 
their statistical theory. On the other hand. Miller (1990) 
2J] notes that for specific initial conditions leading to a 
Gaussian vorticity distribution at statistical equilibrium, 
the mean flow has a linear uJ — tjj relationship similar to 
that obtained from the minimum enstrophy principle. Fi- 
nally, Chavanis & Sommeria (1996) [26| consider a limit 
of strong mixing (or low energy) of the MRS theory in 
which Patp <^ 1 and flnd some connections between max- 
imum entropy states and minimum enstrophy states. In 
that limit, the maximization of entropy at fixed energy, 
circulation and Casimirs becomes equivalent to the mini- 
mization of the macroscopic enstrophy Fj ^' ~ Jzo^dr at 
fixed energy and circulation. This justifies a form of invis- 
cid minimum enstrophy pri nciple in a well-defined limit 
of the statistical theory [5g. This strong mixing limit 
also makes a hierarchy among the Casimir constraints. 
To lowest order in the expansion (Saip ^ 1, leading to 
a linear U — ^ relationship, only the circulation F and 
the microscopic enstrophy T^'^' = J uj-^ dr are important. 
To next orders, leading to nonlinear lj — relationships, 
higher and higher moments F^^f ■ = J w" dr become rele- 
vant. 

Recently, an alternative statistical theory has been pro- 
posed by Ellis, Haven & Turkington (2 002) 1291 and fur- 
ther discussed by Chavanis (2005,2008) [Ifllsil and Cha- 
vanis et al. (2010) [32|. These authors argue that, in real 
situations where the flows are forced and dissipated at 
small scales, the conservation of all the constraints of 
the 2D Euler equation is abusive. They propose to keep 
only the robust constraints (energy and circulation) and 
treat the fragile constraints canonically. This amounts 
to prescribing a prior vorticity distribution instead of the 
Casimirs. This can be viewed as a grand microcanonical 
version of the Miller- Robert-Sommeria (MRS) theory. In 
the Ellis-Haven- Turkington (EHT) approach, the statis- 
tical equilibrium state maximizes a relative entropy (de- 
termined by the prior) while conserving energy and circu- 
lation. The mean flow turns out to maximize a general- 
ized entropy determined by the prior at fixed circulation 
and energy. For a Gaussian prior, the generalized entropy 
is proportional to minus the macroscopic enstrophy. This 
justifies a minimum enstrophy principle from statistical 
mechanics when the constraints are treated canonically 
and the prior is Gaussian. Furthermore, this approach 
allows to go beyond the minimum enstrophy principle by 
considering more complicated priors. 

We can give another interpretation of the EHT ap- 
proach. Indeed, Bouchet (2008) [H notes that the EHT 
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approach provides a sufficient condition of MRS stabil- 
ity. Indeed, it is well-known in statistical mechanics and 
optimization theory [s^l that a solution of a maximiza- 
tion problem is always a solution of a more constrained 
dual maximization problem. Thus, grand microcanonical 
stability (EHT) implies microcanonical stability (MRS). 
Therefore, if the mean flow maximizes a generalized en- 
tropy at fixed circulation and energy, then the corre- 
sponding Gibbs state is a MRS equilibrium state. How- 
ever, the reciprocal is wrong in case of ensemble inequiva- 
lence (between microcanonical and grand microcanonical 
ensembles) that is generic for systems with long-range in- 
teractions nil . The mean flow associated to a MRS equi- 
librium state does not necessarily maximize a generalized 
entropy at fixed circulation and energy. In this sense, the 
minimization of enstrophy at fixed circulation and energy 
provides a sufficient, but not necessary, condition of MRS 
stability. A review of the connections between these dif- 
ferent variational principles has been given recently by 
Chavanis (2009) j3i |. 

In this paper, we shall complement these different ap- 
proaches by considering a sort of intermediate situation 
between all these theories. We argue that, in most physi- 
cal situations, the system is forced and dissipated at small 
scales. In some situations, forcing and dissipation equili- 
brate each other so that the system becomes, in average, 
statistically equivalent to the 2D Euler equation where 
forcing and dissipation are "switched-off" . In particu- 
lar, a quasi stationary state (QSS) can form on a rela- 
tively short timescale. This QSS involves a mean flow 
and fluctuations arount it. We propose to describe this 
state in terms of the statistical mechanics of the 2D Eu- 
ler equation. However, we indirectly take into account 
the effects of forcing and dissipation in the choice of the 
constraints. The energy and the circulation must be ob- 
viously conserved. By contrast, the conservation of all 
the Casimir invariants is abusive and it is likely that 
some Casimir invariants will be destroyed by the forc- 
ing and the dissipation. We argue that some Casimir 
invariants are more relevant than others and that they 
will be selected by the properties of forcing and dissi- 
pation. In that case, the statistical equilibrium state is 
expected to maximize the MRS entropy with these sole 
constraints. It is not our goal here to determine how 
these constraints are selected by the properties of forc- 
ing and dissipation. This is clearly a complicated prob- 
lem that has to be tackled by other methods. We shall 
just use a heuristic approach and consider the situation 
where these relevant constraints are the circulation, the 
energy and the microscopic enstrophy. In other words, 
among all the Casimir constraints, we only consider the 
quadratic one. We do not claim that this is the most 
general situation but simply that it is a case of physi- 
cal interest. We therefore consider the maximization of 
the MRS entropy at fixed circulation, energy and micro- 
scopic enstrophy and provide a detailed study of this (non 
trivial) variational principle. Our maximization principle 
is similar in spirit to the approach of Kraichnan, which 



only considers the conservation of energy and enstrophy. 
However, we remain in physical space and work with the 
MRS entropy for the distribution of vorticity levels while 
Kraichnan works in Fourier space. Furthermor, we take 
into account the conservation of circulation. 

It may be noted that our approach is partly motivated 
by results of experiments carried out in von Karman flows 
[S^ls^l. In fact, the present paper, valid for the 2D Eu- 
ler equations, prepares the ground before considering the 
more complicated (but similar) case of 3D axisymmet- 
ric turbulence treated in [s^. In that later case, we show 
the equivalence between maximum entropy states at fixed 
helicity, angular momentum and microscopic energy and 
minimum macroscopic energy states at fixed helicity and 
angular momentum. Therefore, our two approaches are 
closely related and provide a statistical basis for justify- 
ing the phenomenological minimum enstrophy (2D ) and 
minimum energy (3D axisymmetric) principles [l6|. 

The paper is organized as follows. 

In Sec. |TT1 we discuss some general properties of the 
steady states of the 2D Euler equations. In particular, we 
mention the refined criterion of nonlinear dynamical sta- 
bility given by Ellis et al. p9| based on the maximization 
of a pseudo entropy at fixed circulation and energy. 

In Sec. Illli we recall the phenomenological minimum 
enstrophy principle based on selective viscous decay. 

In Sec. IIV[ we discuss the connection between maxi- 
mum entropy states and minimum enstrophy states. In 
Sec. IIV A[ we consider the maximization of the MRS en- 
tropy at fixed energy, circulation and microscopic enstro- 
phy (energy-enstrophy-circulation statistical mechanics). 
This maximization problem leads to a statistical equilib- 
rium state with Gaussian fluctuations and a mean flow 
characterized by a linear ZJ — ijj relationship. In Sec. 
IIVB| we introduce an equivalent, but simpler, maximiza- 
tion problem based on the maximization of an entropic 
functional of the coarse-grained vorticity at fixed circu- 
lation, energy and microscopic enstrophy. This gener- 
alized entropy has some similarities with the Renyi en- 
tropy (sol. |40|. In Sec. IIV C[ we show the equivalence be- 
tween the maximization of entropy at fixed energy, circu- 
lation and microscopic enstrophy and the minimization of 
macroscopic enstrophy at fixed energy and circulation^^ . 
Therefore, our simplified thermodynamic approach pro- 
vides a justification of the minimum enstrophy principle 
(and Fofonoff flows) from statistical mechanics when only 
the microscopic enstrophy is conserved among the infinite 
class of Casimir constraints. 

We also derive relaxation equations associated with the 
various maximization problems mentioned above (for the 
reader's convenience, the formal derivation of these relax- 
ation equations is postponed to Appendix lD)) . Associated 
with the basic variational problem, we derive a relax- 
ation equation for the vorticity distribution that increases 
the MRS entropy at fixed energy, circulation and micro- 
scopic enstrophy. In that case, the vorticity distribution 
p(r, cr, t) progressively becomes Gaussian and the mean 
flow aJ(r, t) relaxes towards a steady state characterized 
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by a linear oj—iIj relationship. Associated with the simpli- 
fied variational problem, we derive a relaxation equation 
for the mean flow that increases a generalized entropy 
(Renyi-like) while conserving energy, circulation and mi- 
croscopic enstrophy. In that case, the vorticity distribu- 
tion is always Gaussian, even in the out-of-equilibrium 
regime, with a uniform centered variance monotonically 
increasing with time. Associated with the minimum en- 
strophy principle, we derive a relaxation equation for 
the mean flow that dissipates the macroscopic enstrophy 
while conserving energy and circulation. These relax- 
ation equations can serve as numerical algorithms to de- 
termine maximum entropy or minimum enstrophy states 
with appropriate constraints. 

In Sec. |Vl we study minimum enstrophy states at 
fixed energy and circulation in rectangular domains. This 
problem was first considered by Chavanis & Sommeria 
(1996) [2^ who report interesting phase transitions be- 
tween monopoles and dipoles depending on the geometry 
of the domain (e.g., the aspect ratio r of a rectangu- 
lar domain) and on the value of the control parameter 
T'^/E. In particular, for F = 0, they report a transition 
from a monopole to a dipole when the aspect ratio be- 
comes larger than Tc = 1.12. For F 7^ and r < Tc, 
the maximum entropy state is always a monopole and 
for F 7^ and r > Tc, the maximum entropy state 
is a dipole for small T'^/E and a monopole for large 
r'^/E. They also studied the metastability of the solu- 
tions (local entropy maxima) and the possible transition 
between a direct and a reversed monopole. This work has 
been followed recently in different directions. Venaille & 
Bouchet (2009) [4l| have investigated in detail the nature 
of these phase transitions from the viewpoint of statisti- 
cal mechanics. In particular, they showed that the point 
(F = 0,Tc = 1.12) corresponds to a bicritical point sepa- 
rating a microcanonical first order transition line to two 
second order transition lines. Keetels, Clercx & van Hei- 
jst (2009) [i^] have studied minimum enstrophy states 
in rectangular or circular domains with boundary condi- 
tions taking into account the effect of viscosity. Finally, 
Taylor, Borchardt & Helander (2009) [H] have shown 
that the two types of solutions appearing in the study 
of Chavanis & Sommeria [2^ could explain the process 
of "spin-up" discovered by Clercx, Maassen & van Heijst 
(1998) [4^. In Sec. |Vl we recall and complete the main 
results of the approach of Chavanis & Sommeria (26| in 
order to facilitate the discussion of the last section. 

In Sec. IVI[ we use the relaxation equations derived in 
Appendix [D] to illustrate the phase transitions described 
in Sec. |Vl On the basis of the relaxation equations, 
we observe a persistence of unstable states that are sad- 
dle points of entropy. Therefore, we argue that unstable 
saddle points of entropy may play a role in the dynamics 
if the system does not spontaneously generate the per- 
turbations that can destabilize them. We also follow an 
hysteretic cycle as a function of the circulation where 
the hysteresis is due to the robustness of metastable 
states (local entropy maxima). Finally, we briefly de- 



scribe the possibility of transitions between direct and 
reversed monopoles in the presence of stochastic forcing. 

Throughout this paper, we consider the simple case 
of incompressible 2D flows without topography. How- 
ever, the main formalism of the theory can be general- 
ized straightforwardly to account for a topography (or a 
/3-effect) by simply replacing the vorticity lu by the po- 
tential vorticity q = lu + h. Som e ap plications will be 
considered in a companion paper We also assume 

throughout the paper that the domain is of unit area 
A=l. 



II. DYNAMICAL STABILITY OF STEADY 
STATES OF THE 2D EULER EQUATION 

We consider a two-dimensional incompressible and in- 
viscid flow described by the 2D Euler equations 

^ + u-Vw = 0, -AiP = io, (1) 
ot 

where ujz = V x u is the vorticity, ip the stream function 
and u the velocity fleld (z is a unit vector normal to the 
flow). The 2D Euler equation admits an inflnite number 
of steady states of the form 

^ = fW, (2) 

where / is an arbitrary function. They are obtained by 
solving the differential equation 

- AV> = /(^), (3) 

with ip ~ on the domain boundary. 

To determine the dynamical stability of such flows, we 
can make use of the conservation laws of the 2D Euler 
equations. The 2D Euler equations conserve an infinite 
number of integral constraints that are the energy 

J Y'^r^i J cui'dr, (4) 

and the Casiniirs 

Ih = J h{uj) dr, (5) 

where h is an arbitrary function. In particular, all the 
moments of the vorticity F„ = J oj" dr are conserved. 
The first moment F = J w dr is the circulation and the 
second moment F2 = J ui^ dr is the enstrophy. Let us 
consider a special class of Casimirs of the form 

S=- J C{uj) dr, (6) 

where C is a convex function (i.e. C" > 0). These func- 
tionals are called pseudo entropy [1^. Ellis et al. [2^ 
have shown that the maximization problem 

maxjS-M I E[uj] = E, T[u] = F}, (7) 
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determines a steady state of the 2D Euler equation that is 
nonlinearly dynamically stable. This provides a refined 
criterion of nonlinear dynamical stability. The critical 
points of ([7]) are given by the variational principle 

SS - I35E - a5T = 0, (8) 

where /3 and a are Lagrange multipliers. This gives 

C"(a;)==-/3i/;-a^w = F(/3V + a), (9) 

where F{x) = {C'y^{~x). We note that a;'(-i/') = 
— /3/C"(aj) so that is a monotonic function increas- 
ing for /3 < and decreasing for /? > 0. This critical point 
is a steady state of the 2D Euler equation. On the other 
hand, it is a (local) maximum of the pscudo entropy at 
fixed energy and circulation iff 

'^"(^)('^^)^ '^^^\^ j ^'^^'^ 

for all perturbations Slj that conserve energy and circu- 
lation at first order. In that case, it is formally non- 
linearly dynamically stable with respect to the 2D Euler 
equations. This criterion is stronger than the well-known 
Arnol'd theorems that only provide sufficient conditions 
of stability. We note, however, that the refined criterion 
([7]) provides itself just a sufficient condition of nonlinear 
dynamical stability. An even more refined criterion of 
dynamical stability is given by the Kelvin-Arnol'd prin- 
ciple. A review of the connections between these different 
stability criteria has been recently given by Chavanis [s^ . 

III. THE MINIMUM ENSTROPHY PRINCIPLE 

Let us consider the minimization of the enstrophy T2 = 
J up' dr at fixed circulation and energy 

min{r2M |£;[a;] = FM = r}. (11) 

UJ 

The critical points are given by the variational principle 

+ 2(36E + 2aSr = 0, (12) 

where 2/3 and 2a are Lagrange multipliers (the factor 2 
has been introduced for compatibility with the results of 
Sec. WUi . This yields 

w = -At/- = - a- (13) 

This is a steady state of the 2D Euler equation charac- 
terized by a linear uj — ip relationship. On the other hand 
it is a (local) minimum of enstrophy at fixed energy and 
circulation iff 

j {5uj f dr + P j 5u}Si; dr > 0, (14) 

for all perturbations Slj that conserve energy and circu- 
lation at first order. 



There are several interpretations of the minimization 
principle pTjl fsst: 

(i) The minimum enstrophy principle was introduced 
in a phenomenological manner from a selective decay 
principle [l3, [3 [13 • Due to a small viscosity, or other 
source of dissipation, the enstrophy (fragile integral) is 
dissipated while the energy and the circulation (robust 
integrals) are relatively well conserved [s^. It is then ar- 
gued that the system should reach a minimum enstrophy 
state at fixed circulation and energy. Note that there is 
no real justification for this last assumption as discussed 
in [23 ■ The enstrophy could decay without reaching its 
minimum. Furthermore, the minimum potential enstro- 
phy principle is difficult to justify in terms of viscous 
effects for the QG equations (see the Appendix of [45j). 

(ii) If we view F2 as a Casimir of the form ^ , the mini- 
mization principle (jlip is equivalent to the maximization 
principle ^ for the pseudo entropy S' = — iF2. In this 
context, it determines a particular steady state of the 
2D Euler equation that is nonlinearly dynamically stable 
according to the refined stability criterion of Ellis et al. 
[2^. This provides another justification of the minimiza- 
tion problem ([TT]) in relation to the inviscid 2D Euler 
equation. 

(iii) For inviscid flows, the microscopic enstrophy 
1^2" ~ J Ld^ dr is conserved by the 2D Euler equation 
but the macroscopic enstrophy F2^' = JZJ^dr calculated 
with the coarse-grained vorticity decreases as enstrophy 
is lost in the fluctuations. Indeed, by Schwartz inequal- 
ity: T^'^' = JiJ^dr< J uj'^ dr = Fj'^'. By contrast, the 
energy E ~ ^ J ZJ'ip dr and the circulation F = JuJdr 
calculated with the coarse-grained vorticity are approx- 
imately conserved. This suggests an inviscid minimum 
enstrophy principle based on the minimization of macro- 
scopic enstrophy at flxed energy and circulation [26l |. In 
this case, selective decay is due to the operation of coarse- 
graining, not viscosity. 

In the following, we shall discuss some connections be- 
tween the minimum enstrophy principle and the maxi- 
mum entropy principle. 

IV. CONNECTION BETWEEN MAXIMUM 
ENTROPY STATES AND MINIMUM 
ENSTROPHY STATES 

A. Energy-enstrophy-circulation statistical theory 

Starting from a generically unstable or unsteady initial 
condition, the 2D Euler equations are known to develop a 
complicated mixing process leading ultimately to a quasi 
stationary state (QSS), a vortex or a jet, on the coarse- 
grained scale. In order to describe this QSS and the 
fluctuations around it, we must introduce a probabilis- 
tic description. Let us introduce the density probability 
p(r, a) of flnding the vorticity level a; = cr at position r. 
Then, the local moments of vorticity are w" = / per" dr. 
In the statistical mechanics approach of Miller-Robert- 
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Sommcria (2J, |25| , assuming that the system is strictly 
described by the 2D Euler equation (no forcing and no 
dissipation), the statistical equilibrium state is expected 
to maximize the mixing entropy 



where a;2 



S[p] 



pin p drda, 



while conserving all the invariants (energy and Casimirs) 
of the 2D Euler equation. This forms the standard MRS 
theory. 

In the case of flows that arc forced and dissipated at 
small scales, one may argue that forcing and dissipation 
will compensate each other in average so that the system 
will again achieve a QSS that is a stationary solution of 
the 2D Euler equation. This QSS will be selected by 
forcing and dissi pat ion. This fact is vindicated both in 
experiments [36l. |37| and numerical simulations [i^. In 
order to describe the fluctuations around this state, one 
needs to go one step further and obtain the vorticity dis- 
tribution. An idea is to keep the framework of the sta- 
tistical theory but argue that forcing and dissipation will 
alter the constraints. More precisely, forcing and dis- 
sipation will select some particular relevant constraints 
among all the invariants of the ideal 2D Euler equation. 
These constraints will determine the mean flow and the 
fluctuations around it. For example, we argue that there 
exists physical situations in which only the conservation 
of energy, circulation and microscopic enstrophy are rele- 
vant (and of course the normalization condition). We do 
not claim that such situations are universal but simply 
that they happen in some cases of physical interest. This 
seems to be the case for example in some oceanic situ- 
ations [6l-[20j (see also 47 1 for recent studies). We shall 
therefore consider the maximization problem 

max{5[p] I E, r, tI-'-, J pda = 1}, (16) 



where 



i? = - / uj^jj dr ^ - I i/jpa drda, 



T = I Lodv ~ / p<7 drda, 



j.-dr = j 



pa" drda. 



(17) 



(18) 



(19) 



The last constraint ([T9)) will be called the microscopic (or 
fine-grained) enstrophy because it takes into account the 
fluctuations of the vorticity w. It is different from the 
macroscopic (or coarse-grained) enstrophy 



^ 2 



uP' dr. 



which ignores these fluctuations. We have 



(20) 



(21) 



uj is the local centered variance of 

fluct 



the vorticity. The fluctuations of enstrophy are 
J ijj2dr. In our terminology, the enstrophy will be called a 
fragile constraint because it cannot be expressed in terms 
of the coarse-grained field since a;^ ^ Zp' . While the 
(15) microscopic enstrophy V^^' is conserved, the macroscopic 
enstrophy ^' is not conserved and decays. By contrast, 



the energy p?|) and the circulation (fT8)) will be called 
robust constraints because they can be expressed in terms 
of the coarse-grained fields. We shall come back to this 
important distinction in Sec. IIV CI 

The maximization problem (|16p is what we shall call 
"the energy-enstrophy-circulation statistical theory" , or 
simply, "the statistical theory" in this paper. We note 
that a solution of (fT6)) is always a MRS statistical equi- 
librium state, but the reciprocal is wrong in case of en- 
semble inequivalence because we have relaxed some con- 
straints (the other Casimirs). Here, we keep the robust 
constraints E and F and only one fragile constraint Fj'^', 
the quadratic one. We assume that these constraints are 
selected by the properties of forcing and dissipation. As 
we shall see, this assumption leads to Gaussian fluctua- 
tions and a mean flow characterized by a linear w — tfj 
relationship. In principle, we can obtain more complex 
fluctuations and more complex mean flows (characterized 
by nonlinear uj — relationships) by keeping more and 
more flne-grained moments F^'^j^ among the constraints. 
This can be a practical way to go beyond the Gaussian 
approximation. However, the Gaussian approximation, 
leading to a linear uJ—tp relationship, is already an inter- 
esting problem presenting rich bifurcations (because the 
energy constraint is nonlinear) (26j . so that we shall stick 
to that situation. 

The critical points of (|16p are solution of the varia- 
tional principle 

SS ~ /3SE - aSr - a2<5F{-^- - J C{r)6 (^J pda^ dr = 0, 



(22) 



where /3, a, a2 and C(r) are Lagrange multipliers. This 
yields the Gibbs state 



P(r,cr) 



Z{r) 



(23) 



where the "partition function" is determined via the nor- 
malization condition 



Z(r) 



da. 



(24) 



Therefore, in this approach, the distribution p(r, a) of 
the fluctuations of vorticity is Gaussian and the centered 
variance of the vorticity W2(r) is uniform 



1 



UJ2 — UJ"" — UJ = — il2. 

2a2 



(25) 



On the other hand, the mean flow is given by 

uj = -il2{f3ij + a). (26) 
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This is a steady state of the 2D Euler equation eharae- 
terized by a hnear 115 — ijj relationship. Then, the Gibbs 
state can be rewritten 



Since cj2 
get 



P(i',cr) 



1 



(27) 



^/2^dh 

is uniform at statistical equilibrium, we 



Since the specification of a;(r) determines E and F, this 
is equivalent to maximizing S at fixed r^'^-, J pda = I 
and w(r) = J pa da. Writing the variational problem as 



5S 



X{r)5 pada \ dr 



J C(r)(5 (^J pda^ dr = 0, 



(31) 



(28) 



we obtain 



Finally, a critical point of 
at fixed E, T and T^'^' iff 



Ipj) is an entropy maximum 



^^drda^lp 
P 2^ 



(5aJ(5V dr <0 (29) 



for all perturbations 6p that conserve energy, circulation, 
microscopic enstrophy and normalization at first order 
(the proof is similar to the one given in [s^ for related 
maximization problems). 

Remark: From Eqs. ^15\\ and (|27)) . we easily find that 
S = ilnf22. Then, using Eq. (1^ . we conclude that, at 
equilibrium, the entropy is given by 



S 



1 



■In 



(30) 



Therefore, if there exists several local entropy maxima 
(metastable states) for the same values of the constraints 



E, F and F 



f-g- 



with the smallest enstrophy T^'^ 



the maximum entropy state is the one 
This is a first re- 
sult showing connections between maximum entropy and 
minimum enstrophy principles. However, this equilib- 
rium result does not prove that the maximization of the 
entropy functional S[p] at fixed E, F and ^2^' is equiv- 
alent to the minimization of the enstrophy functional 
F2^' [ZU] at fixed E and F (e.g. an entropy maximum 
could be a saddle point of enstrophy) . This equivalence 
will be shown in See. IIV CI 

Remark: the maximization problem (|16p also arises in 
the study of Kazantsev et al. when the relaxation 
equations associated with the MRS statistical the ory are 
closed by using a Gaussian approximation (see [32] for 
details). 



Pi(r,cr) 



with 



2a2 



(32) 



(33) 



Note that the centered variance of the vorticity lo2{v) = 
U,2 is uniform. Equation (j33p also implies that a2 
must be positive. We check that pi is a global entropy 
maximum with the previous constraints since 5^S = 

— J -^^p drda < (the constraints are linear in p so their 
second variations vanish) . Using the optimal distribution 
(|32p . we can express the entropy (|15|1 in terms of w writ- 
ing S\lJ] = S[pi]. After straightforward calculations we 
obtain 



5= iln02, 



(34) 



up to some unimportant constant terms. Note that f22 is 
determined by the constraint on the microscopic enstro- 
phy which leads to 



n2 = H'- 



dr. 



(35) 



The second term is the macroscopic enstrophy associated 



with the coarse-grained flow: 



pc.g. 
^ 2 



/ uj^ dr. The rela- 



tion psp can be used to express the entropy (|34|) in terms 
of ZJ alone. 

(ii) Second step: we now have to solve the maximiza- 
tion problem 



max{S[uJ]\E, F, F^'^''}, 



(36) 



B. An equivalent but simpler variational principle 



with 



The maximization problem (|16p is difficult to solve, 
especially regarding the stability condition (f29| , because 
we have to deal with a distribution p(r,a). We shall 
introduce here an equivalent but simpler maximization 
problem by "projecting" the distribution on a smaller 
subspace. To solve the maximization problem (|16p . we 
can proceed in two steps [60j: 

(i) First step: We first maximize S at fixed E, 
J pda — 1 and a given vorticity profile SJ(r) — J pa da. 



F, F^-^- 



■In 



E = 



rf-g- 

^ 2 



dr 



Ljip dr, 



F = / uJdr. 



(37) 



(38) 



(39) 
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The functional (|37|) might be called a generalized entropy. 
Interestingly, it resembles the Renyi entropy [s^ \^ . 

(in) Conclusion: finally, the solution of (|16|) is given by 
Eq. (1321) where oJ is determined by Therefore, 
and (p6| are equivalent but ((36)) is easier to solve because 
it is expressed in terms of uj(r) while p6p is expressed in 
terms of p(r, a). 

Up to second order, the variations of entropy p7p are 
given by 



AS*: 



1 



LoSuj dr - 



[SuY dr 



jSu! dr 



(40) 

where V,2 is given by Eq. (|35l). Considering the first 
order variations of the entropy, the critical points of p6p 
are determined by the variational problem 



This yields 



6S - PSE - a5T = 0. 



w = -172(/3V + a). 



(41) 



(42) 



We recover Eq. ([26)) for the mean flow. Combined with 
Eq. p2p . we recover the Gibbs state ([231) ■ Considering 
now the second order variations of the entropy ([M)) , we 
find that a critical point of p6p is a maximum of entropy 
at fixed energy, circulation and microscopic enstrophy iff 

1 



2^2 I ^^^^^'^^^^ ^ Si^Si^dr 



1 



uiSu) dr ] < 0, 



(43) 



for all perturbations SZJ that conserve circulation and en- 
ergy at first order (the conservation of microscopic en- 
strophy is automatically taken into account in our for- 
mulation). This stability condition is equivalent to Eq. 
()29p but much simpler because it depends only on the 
perturbation SoJ instead of the perturbation of the full 
distribution 6p. In fact, the stability condition ((43| can 
be simplified further. Indeed, using Eq. ()42p . we find 
that the term in parenthesis can be written 



LoSu) dr ~ — fl2 / {P'>p + a)duj dr 



(44) 



and it vanishes since the energy and the circulation arc 
conserved at first order so that SE = / i/jSlj dr ~ and 
ST = J Suj dr = 0. Therefore, a critical point of ((36)) is 
a maximum of entropy at fixed energy, circulation and 
microscopic enstrophy iff 



(45) 



" 2^ / ^^^^^ " f / ^^^^'^'^ °' 

for all perturbations SZJ that conserve circulation and en- 
ergy at first order. In fact, this stability condition can 
be obtained more rapidly if we remark that the maxi- 
mization problem p6p is equivalent to the minimization 
of the macroscopic enstrophy at fixed energy and circu- 
lation (see Sec. lIVCp . 



C. Equivalence with the minimum enstrophy 
principle 

Since ln(x) is a monotonically increasing function, it 
is clear that the maximization problem p6p is equivalent 
to 



with 



max{S'[cU] \E,r}, 



= / dr. 



- I ujiPdr, 



r = / Udr. 



(46) 



(47) 



(48) 



(49) 



(50) 



The functional S of the coarse-grained vorticity ZJ is 
called a "generalized entropy" . It is proportional to the 
opposite of the coarse-grained enstrophy. We have the 
equivalences 



(51) 



Therefore, the maximization of MRS entropy at fixed en- 
ergy, circulation and microscopic enstrophy is equivalent 
to the minimization of macroscopic enstrophy at fixed 
energy and circulation. The solution of ()16p is given by 
Eq. ((32)) where ZU is determined by ((46)) and ^2 by Eq. 
((55)) . Therefore, ()16p and ()46p are equivalent but ((^5)) is 
easier to solve because it is expressed in terms of ZU while 
(|16p is expressed in terms of p. This provides a justifica- 
tion of the coarse-grained minimum enstrophy principle 
in terms of statistical mechanics when only the micro- 
scopic enstrophy is conserved among the Casimirs. Note 
that, according to ((7)), the principle ((46)) also assures that 
the mean fiow associated with the statistical equilibrium 
state ()16p is nonlinearly dynamically stable with respect 
to the 2D Euler equation. 

The critical points of ((46)) are given by the variational 
problem 



This yields 



5S - PSE - a5T = 0. 



uj ~ — ~ Q^- 



(52) 



(53) 



This returns Eq. ((26)) for the mean fiow (up to a trivial 
redefinition of (3 and a). Together with Eq. ((32)), this 
returns the Gibbs state ((23)) . On the other hand, this 
state is a maximum of S at fixed E and F iff 

- i / {Suj)^ * " f / ^^^'^ '^^ < (5"^) 
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for all perturbations StJ that conserve circulation and en- 
ergy at first order. This is equivalent to the criterion ()45|) 
as it should. 

We have thus shown the equivalence between the maxi- 
mization of MRS entropy at fixed energy, circulation and 
fine-grained enstrophy with the minimization of coarse- 
grained enstrophy at fixed energy and circulation. This 
equivalence has been shown here for global maximiza- 
tion. In Appendix |X1 we prove the equivalence for local 
maximization by showing that the stability criteria ()29|) 
and (|5^ are equivalent. 



D. Equivalence with a grand microcanonical 
ensemble 



In the basic maximization problem (jl6|) . the fine- 
grained enstrophy is treated as a constraint. Let us in- 
troduce a grand microcanonical ensemble by making a 
Legendre transform of the entropy with respect to this 
fragile constraint T^'^' [s^. We thus introduce the func- 



tional So 



and the maximization problem 



max{S'g[p] \ E, T, / pda = 1}. 
p J 



(55) 



A solution of ((55|) is always a solution of the more con- 
strained dual problem (|16p but the reciprocal is wrong in 
case of "ensemble inequivalence" [Hj. In the present case, 
however, we shall show that the microcanonical ensem- 
ble ([T6| and the grand microcanonical ensemble ([SS]) are 
equivalent. This is because only a quadratic constraint 
(enstrophy) is involved. 

To solve the maximization problem (j55p we can pro- 
ceed in two steps. We first maximize Sg at fixed E, F, 
J pda — 1 and ZU(r) = J pa da. This is equivalent to 
maximizing Sg at fixed J pda = 1 and ^(r) = J pa da, 
and this leads to the optimal distribution (|32p where 
^2 = 1/(2q;2) is now fixed. This is clearly the global 
maximum of Sg with the previous constraints. Using 
this optimal distribution, we can now express the func- 
tional Sg in terms of oJ by writing S\uJ] = Sg[pi]. After 
straightforward calculations, we obtain 



(|46p . Since we have proven previously that (|46p is equiv- 
alent to the microcanonical variational principle (|16p , we 
conclude that (fT6|) and (|55|) are equivalent. 

Remark: the grand microcanonical ensemble (|55p cor- 
responds to the EHT approach with a Gaussian prior 

iilllil. 



E. Connection between different variational 
principles 

Let us finally discuss the relationship between our ap- 
proach, Naso-Chavanis-Dubrulle (NCD), and the ones 
proposed by Miller-Robert-Sommeria (MRS) and Ellis- 
Haven- Turkington (EHT). To that purpose, we shall 
make the connection between the variational principles 



(MRS) : max{S[p] \ E, F, f{|i, y pda = 1}, (58) 
(EHT) : inax{S^ [p] \E,T, j pda = 1}, (59) 
(NCD) : miix{S[p]\E,T,rl-^\ J pda ^ 1}. (60) 



(MaxS) : max{S[ZJ]\E,T}. 

(MinFa) : niin{F^-» [cJ] | F}, 
where the functionals are 

S[p\ = — J p{r, a) In p{r, a) drda, 



^xip] = - / P(i",cr)ln 



p(r,o-) 



drda. 



(61) 
(62) 

(63) 
(64) 



Sa 



2a 



pc.g. 
2 



(56) 



up to some constant terms (recall that is a fixed pa- 
rameter in the present situation). In the second step, we 
have to solve the maximization problem 



max{^[ZU] I E, F}. 



Finally, the solution of ([55)1 is given by Eq. ([5^ where 
U is determined by (j57p . Therefore, the variational prin- 
ciple ([55]) is equivalent to (|57)) . That this is true also 



for local maximization is shown in Appendix B of [32 1 
(in a more general situation). On the other hand, since 
fl2 > 0, the maximization problem (|57p is equivalent to 



S[u 



C{iO) dr, 



pc.g. 
^ 2 



' dr. 



(65) 



(66) 



(57) with x(o-) 



exp(-X;„>iQ;„o-") and C(a;) = 

-/"[(lnx)']"^(-a;)dx where x($) = J xia)e-'"^ da. 

In the framework of the MRS approach where all the 
Casimirs are conserved, the maximization of a "gener- 
alized entropy" S\lJ] at fixed energy and circulation pro- 
vides a su ffici ent condition of MRS thermodynamical sta- 
bility [3^ l33l Issj . However, the reciprocal is wrong in 
case of "ensemble inequivalence" between microcanonical 
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and grand microcanonical ensembles. Indeed, the coarse- 
grained vorticity field ZU(r) associated with a MRS ther- 
modynamical equilibrium (i.e. a maximum of entropy 
S[p\ at fixed energy, circulation and Casimirs) is not nec- 
essarily a maximum of generalized entropy S[uj] at fixed 
energy and circulation (it can be a saddle point of gen- 
eralized entropy at fixed energy and circulation). 

In the framework of the EHT approach where the con- 
servation of the Casimirs is replaced by the specifica- 
tion of a prior vorticity distribution (i.e. the Casimirs 
are treated canonically) , the maximization of a general- 
ized entropy S\Uj] at fixed energy and circulation provides 
a necessary and sufficient condition of EHT thermody- 
namical stability [l^, [HI HI]- Indeed, a vorticity dis- 
tribution p{v, a) is a EHT thermodynamical equilibrium 
(i.e. a maximum of relative entropy Sy\p] at fixed energy 
and circulation) if and only if the corresponding coarse- 
grained vorticity field w(r) is a maximum of generalized 
entropy S\ijj] at fixed energy and circulation. 

Thus, we symbolically have 

(MRS) (EHT) ^ (MaxS) (67) 

Let us now specialize on the case of Gaussian distribu- 
tions. 

In the framework of the MRS approach where all the 
Casimirs are conserved, the minimization of macroscopic 
enstrophy ^' [w] at fixed energy and circulation pro- 
vides a sufficient condition of MRS thermodynamical sta- 
bility for initial conditions leading to a Gaussian vorticity 
distribution at equilibrium [3^ Issl. Issj . However, the re- 
ciprocal is wrong in case of "ensemble inequivalence" , i.e. 
the coarse-grained vorticity field uJ(r) associated with a 
MRS thermodynamical equilibrium with Gaussian vor- 
ticity distribution is not necessarily a minimum enstro- 
phy state (it can be a saddle point of macroscopic enstro- 
phy at fixed energy and circulation). 

In the framework of the EHT approach where the con- 
servation of the Casimirs is replaced by the specifica- 
tion of a prior vorticity distribution, the minimization 
of macroscopic enstrophy r2^ [w] at fixed energy and 
circulation provides a necessary and sufficient condition 
of EHT thermodynamical stability for a Gaussian prior 
[Hill Hi, i.e. a vorticity distribution p(r, cr) is a EHT 
thermodynamical equilibrium with a Gaussian prior if 
and only if the corresponding coarse-grained vorticity 
field il;(r) is a minimum of macroscopic enstrophy ^' \Uj\ 
at fixed energy and circulation. 

In the framework of the NCD approach where only 
the microscopic enstrophy F2'®' [p] is conserved among the 
Casimir constraints, the minimization of macroscopic en- 
strophy Fj ^'[oJ] at fixed energy and circulation provides a 
necessary and sufficient condition of NCD thermodynam- 
ical stability. Indeed, a vorticity distribution p{v, cr) is a 
NCD thermodynamical equilibrium (i.e. a maximum of 
entropy S[p\ at fixed energy, circulation and microscopic 
enstrophy) if and only if the corresponding coarse-grained 
vorticity field aJ(r) is a minimum of macroscopic enstro- 
phy T2^' [aJ] at fixed energy and circulation. 



Thus, we symbolically have 

(MRS) ^ (EHT) ^ (MinFs) ^ (NCD) (68) 

Remark 1: the EHT and NCD approaches provide suf- 
ficient conditions of MRS stability. They are valuable 
in that respect as they are simpler to solve. They may 
also have a deeper physical meaning as discussed in the 
introduction. 

Remark 2: the equivalence between (EHT) for a Gaus- 
sian prior and (NCD) is essentially coincidental because 
these variational problems are physically very different. 
In particular, this agreement is only valid for a Gaussian 
prior and does not extend to more general cases. 

Finally, for completeness, we mention similar results 
obtained in [ssj for axisymmetric flows. In the frame- 
work of the Naso-Monchaux-Chavanis-DubruUe (NMCD) 
approach where only the microscopic energy E^'^' [p\ is 
conserved among the Casimir constraints, the minimiza- 
tion of macroscopic energy £"^f • \a] at fixed hclicity and 
angular momentum provides a necessary and sufficient 
condition of NMCD thermodynamical stability, i.e. a 
distribution of angular momentum p[t, rf) is a maximum 
of entropy S[p] at fixed helicity, angular momentum and 
microscopic energy E^'^' [p] if and only if the correspond- 
ing coarse-grained angular momentum distribution ^(r) 
is a minimum of macroscopic energy iJ'^f - [a] at fixed he- 
licity and angular momentum. 

Thus, we symbolically have 

(NMCD) ^ (Mini;) (69) 



V. PHASE TRANSITIONS IN 2D EULER 
FLOWS 

A. Minimum enstrophy states 

In this section we study the maximization problem 

m3.-K{Sp] \E,T}, (70) 

UJ 

where S = — ^ / dr is the neg-enstrophy (the oppo- 
site of the enstrophy), E ~ ^Joj^jjdr the energy and 
F = J u} dr the circulation. The maximization problem 
(|70p can be interpreted as: (i) a criterion of nonlinear 
dynamical stability with respect to the 2D Euler equa- 
tion (Sec. nil, (ii) a phenomenological minimum enstro- 
phy principle (Sec. IIIip . (iii) a sufficient condition of 
MRS thermodynamical stability [H, HH, (iv) a neces- 
sary and sufficient condition of EHT thermodynamical 
stability for a Gaussian prior [sil, [13], (v) a necessary 
and sufficient condition of thermodynamical stability in 
the energy-enstrophy-circulation statistical theory where 
only the microscopic enstrophy is conserved among the 
Casimirs (Sec. IIV[) . For simplicity and convenience, we 
shall call S the entropy. 
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We write the variational principle for the first order 
variations as 



SS ~ P5E - a6T = 0, 



(71) 



where /3 and a are Lagrange multipliers. This yields a 
linear uj — ip relationship 



cj = —Alp = —/Sip — a. 



(72) 



As before, we assume that the area of the domain is unity 
and we set {X) = J X dr. Taking the space average of 
Eq. ([72|l. we obtain T = — /3('0) — a so that the foregoing 
equation can be rewritten 



- AV' + /?V- = r + /j(V'), 



(73) 



with ■(/) = on the domain boundary [6l|. This is the 
fundamental differential equation of the problem. The 
energy and the entropy can then be expressed as 



S 



(74) 



(75) 



We shall study the maximization problem (1701) bv adapt- 
ing the approach of Chavanis & Sommeria [26| to this 
specific situation (these authors studied a related but 
not exactly equivalent problem). We will see that the 
structure of the problem depends on a unique control 
parameter (26j : 



A 



'2E 



(76) 



We note that the maximization problem ([701) has been 
studied recently by Venaille & Bouchet [41| by using a 
different theoretical treatment. They performed a de- 
tailed analysis of the phase transitions associated with 
([70]) in the context of statistical mechanics, emphasiz- 
ing in particular the notion of ensemble inequivalence. 
However, their approach is very abstract. Our study is 
more direct and can offer a complementary discussion of 
the problem. The maximization problem (j70p has also 
been studied recently by Keetels et al. [l^l with different 
boundary conditions adapted to viscous flows. 



B. The bifurcation diagram 

In this section, we apply the methodology developed 
by Chavanis & Sommeria [2Q| . This methodology is rela- 
tively general: it is valid for an arbitrary domain and for 
an arbitrary linear operator. However, for illustration, 
we shall consider the Laplacian operator and a rectangu- 
lar domain. 



1. The eigenmodes 
We first assume that 



(77) 



corresponding to a = 0. In that case, the differential 
equation ([75)) becomes 



At/; + ,3V' = 0, 



(78) 



with V' = on the domain boundary. Using the results 
of Appendix [HI Eq. ([78[) has solutions only for /3 = /3„i„ 
(eigenvalues) and the corresponding solutions (eigenfunc- 
tions) are 



V' = 



2E 

~ Pnii 



1/2 



(79) 



where we have used the energy constraint ()74[) to deter- 
mine the normalization constant. Substituting this result 
in Eq. (|77p . we find that these solutions exist only for 
A = Amn with 



A' 



n 



(80) 



For the eigenmodes {ipmn) — (m or n even), we find 
A = and for the eigenmodes {tpmn) 7^ (m and n odd). 



we find A^ = A". 



r? 



^0. 



2. The solutions of the continuum 



We now assume that F + /3 (t/j) ^ and we define 



(81) 



" F + l3{xP) 
In that case, the differential equation ([73[) becomes 

-At/) + /?(/) = 1, (82) 

with = on the domain boundary. We also assume 
that (3 7^ Pmn- In that case, Eq. ([5^ has a unique 
solution that can be obtained by expanding on the 
eigenmodes. We get 



E('0?riTi) , 



(83) 



where only the modes with {ipmn) 7^ are "excited" . 

For F 7^ 0, taking the average of Eq. ()81|) and solving 
for {tp), we obtain {tp) = T{(p)/{1- P3{(p)). Therefore, the 
solution of Eq. ([73[) is 



T4> 



1 - ■ 



(84) 



Substituting this solution in the energy constraint (|74p . 
we obtain the "equation of state" : 



(85) 
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This equation determines /3 as a function of A. In par- 
ticular, it determines the caloric curve /3{E) for a given 
value of r 7^ 0. Note that the equation of state involves 
the important function [26| : 

F{I3) = /3(0) - 1. (86) 

For r = 0, the solution of Eq. is 

V' = /SW^. (87) 

Taking the space average of this relation, we find that 
this solution exists only for a discrete set of temperatures 
/3 = pi""^ satisfying F{pi''^) = 0. We shall note = pi^'> 
the largest of these solutions. Substituting Eq. ([57]) in 
the energy constraint ([74)) . wc find that the amplitude 
(■0) is determined by 

S = -^/?'W'((0'> (88) 

Of course, the case A = is also a limit case of the 
equation of state 

3. The mixed solutions 
For /? — > /3„in with m and n odd ((V'mn) 7^ 0), we find 

from Eq. dSS]) that </> ^ ('0mn)V'mn/(/3 - /3mn) -5> +0O 

leading to A AJ^„ and ip (2£'//3m„)^/^Vrnri- There- 
fore, we recover the eigenfunction '(/'mn as a limit case. 
The eigenfunctions with non vanishing average value are 
therefore contained on the continuum branch. 

For /? = with m or n even ((V'mn) = 0), the 

solution of Eq. ((5^ is not unique. It corresponds to the 
mixed solutions 

<^ = ft o ^ V'm'n' + XmnV'mn, (89) 

m'n' 

where Xmn is determined by the energy constraint (more 
precisely, it can be related to A by substituting Eq. ([M)) 
in Eq. (|85p where now (3 = Pmn)- These solutions form a 
plateau at fixed temperature /3 = /3mn- For Xmn — >■ +00, 
we recover the pure eigenmode i^mn that exists at A = 
and for Xmn = 0, we connect the branch of continuum 
solutions at A = A^„. 

The general bifurcation diagram showing the eigen- 
modes, the solutions of the continuum and the mixed 
solutions is represented in Fig. 2 of [1^ (see also Figs. [1] 
and [2] below). 

C. The geometry induced monopole/dipole 
transition 

For a given value of the control parameter A, we can 
have an infinite number of solutions to Eq. (|73p |2q | . We 
can now use the entropy ([75|) to select the most probable 
state (maximum entropy state) among all these solutions. 



For the eigenmodes, the entropy takes the simple form 
S/E = Pmn- In particular, for the eigenmodes '0„m with 
m or n even {{ipmn) = 0) that exist only for A = 0, we 
have 

^(A = 0) (90) 

For a rectangular domain elongated in the x direction, 
the eigenmode with the highest entropy at F = is the 
dipole {m,n) = (2,1) with temperature /32i(''")- There- 
fore, the maximum entropy state (or the minimum en- 
strophy state) corresponds to the mode with the largest 
scale. The modes with smaller scales (to, n large) have 
lower entropy (higher enstrophy). Therefore, the max- 
imum entropy and minimum enstrophy principles se- 
lect the large-scale structures among the infinite class 
of steady states of the 2D Euler equation. This is a man- 
ifestation of the inverse cascade process. 

For the solutions of the continuum, the entropy can be 
written 

For A = 0, this expression reduces to 

|(A = 0) = M'^ (92) 

The solution with highest entropy is the monopole with 
temperature /3*(r). 

For A = 0; we have a competition between the 
monopole /3*(r) (continuum branch) and the dipole 
P2i{t) (eigenmode) [62| . We must therefore compare their 
entropy (or equivalently their inverse temperature) to se- 
lect the maximum entropy state. As shown in Chavanis 
& Sommeria [2^, this selection depends on the geome- 
try of the domain. In a rectangular domain, it is found 
that the monopole has the highest entropy (/3* > /32i) 
for T < Tc = 1.12 while the dipole dominates (/32i > /?*) 
for T > Tc = 1.12. More generally, it can be shown that 
the entropy is a monotonically increasing function of the 
inverse temperature (for a fixed value of A). Therefore, 
at any A, the maximum entrop y s tate is the one with 
the highest inverse temperature |26| . The series of equi- 
libria /3(A) is represented in Figs. [T] and [5] for a square 
domain and for a rectangular domain of aspect ratio 2, 
respectively. For r < Tc, the maximum entropy state is 
the direct monopole for any value of A. For r > Tc, the 
maximum entropy state is the dipole for A^ < (Aj^)^ and 
the direct monopole for A^ > (Aji)^. 

D. Stability analysis and ensemble inequivalence 

For given values of E and F, there can exist different 
critical points of entropy S (canceling its first order vari- 
ations). They are solutions of the differential equation 
([75)1 . For sufficiently small A, there exists an infinity of 
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FIG. 1: Series of equilibria in a square domain (t = 1 < Tc). 
In that case max{/3«, /32i} = P»- The maximum entropy state 
is the direct monopole (for F > the vorticity is positive at 
the center and negative at the periphery (MP); for F < this 
is the opposite (MN)) for any value of For < {A'2if, 
the reversed monopole is metastable (local entropy maximum) 
as discussed in Sec. IV Dl Note that the metastable states 
have negative specific heats C = |f = ^^£;^^%^ < 0. 
For F = 0, the direct and reversed monopoles have the same 
entropy. For fixed F, the caloric curve 13(E) does not present 
any phase transition (see Sec. IV Fp . The vorticity profiles 
are plotted for F > 0. The red colour corresponds to positive 
values of the vorticity and the blue colour to negative values. 




FIG. 2: Series of equilibria in a rectangular domain with as- 
pect ratio r = 2 > Tc. In that case max{/3t, /32i} = /32i. The 
maximum entropy state is the dipole for A^ < (A2i)'^ and 
the direct monopole for A^ > (A2i)^ (the reversed monopoles 
are unstable). For F 7^ 0, the caloric curve l3{E) presents a 
second order phase transition marked by the discontinuity of 
§§{E) at E = E'21 as discussed in Sec. IVF] 



solutions [26|. In the last section, we have compared the 
value of the entropy of these different solutions in order to 
select the maximum entropy state. However, a more pre- 
cise study should determine which solutions correspond 
to global entropy maxima, local entropy maxima and sad- 
dle points. Saddle points of entropy are unstable and 
should be rejected in principle (see, however Sec. IVIBI) . 



By contrast, local entropy maxima (metastable states) 
can be long-lived for systems with long-range interac- 
tions. In practice, they are as much relevant as global 
entropy maxima (stable states). In the following, using 
an approach very close to the one followed by Chavanis 
& Sommeria |2al (but not exactly equivalent since the 
variational problems differ), we determine sufficient con- 
ditions of instability. This will eliminate a large class of 
solutions that are unstable saddle points of entropy and 
give the form of the perturbations that destabilize them. 
The remaining solutions are either stable or metastable. 

A critical point of entropy at fixed energy and circula- 
tion is a (local) maximum iff 



5' J 



5uj5tl)dr < 0, 



(93) 



for all perturbations that conserve energy and circulation 
at first order: 5E = 5T ^ 0. 

(i) We first show that all the solutions with (3 < f32i 
arc unstable (saddle points). To that purpose, we con- 
sider a perturbation of the form Suj = ■(/'2i(r)- The corre- 
sponding stream function is Sip = — ^V'2i(r). For this 
perturbation, it is clear that JT = J Sui dr = since 
("021) = 0. Furthermore, 5E ~ J il'Sio dr = since tp2i is 
orthogonal to the other eigenmodes and to the solutions 
of the continuum (as they involve a summation (|83p on 
the eigenmodes with non zero average that are orthogo- 
nal to '021 )• Finally, a simple calculation shows that 



6' J 



/321 



1 > 0. 



(94) 



We have thus found a particular perturbation that in- 
creases the entropy at fixed energy and circulation. 
Therefore, the states with /3 < /32i are unstable (saddle 
points). 

(ii) We now show that if /32i < /?* the mode 1^21 exist- 
ing at A = is unstable (saddle point). To that purpose, 
we consider a perturbation of the form Suj = 1 — /3*c/)* 
(where is the solution of Eq. (|82p corresponding to 
/? = /?*). The corresponding stream function is Sip = (f>^. 
For this perturbation, it is clear that ST = J Suj dr = 
since 1 — /3*(0*) = 0. Furthermore, SE = J ip2iSujdr = 
since (V'21) = and ^p2l is perpendicular to as ex- 
plained previously. Finally, after some simple algebra 
using 1 — /3*(0*) — 0, we get 

S^J = liP* - /32i)((0,) - PM)) > 0, (95) 

(the last term in parenthesis is positive as shown in Ap- 
pendix|B]). We have thus found a particular perturbation 
that increases the entropy at fixed energy and circulation. 
Therefore, if /32i < /3* the mode -021 existing at A = 
is unstable. By continuity, the mixed solutions forming 
a plateau at /3 = /32i are also unstable if P21 < P* since 
the two ends of the plateau are unstable. 

The maximization problem (|70|) corresponds to a con- 
dition of microcanonical stability which is relevant to our 
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problem since the circulation and the energy are con- 
served by the 2D Euler equation. However, it can be 
convenient to establish criteria of canonical and grand 
canonical stability. Indeed, the solution of a maximiza- 
tion problem is always solution of a more constrained 
dual maximization problem, but the reciprocal is wrong 
in case of ensemble inequivalence that is generic for sys- 
tems with long-range interactions 3J]. Therefore, con- 
ditions of canonical and grand canonical stability pro- 
vide only sujficient conditions of microcanonical stabil- 
ity: grand canonical stability implies canonical stabil- 
ity which itself implies microcanonical stability (see, e.g., 
[35l|). This problem of ensemble inequivalence has been 
studied in detail by Venaille & Bouchet [4l[ and we briefly 
discuss it again bringing some complements regarding the 
metastable states (that arc not considered in 
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Considering the grand canonical ensemble, we have to 
maximize the grand potential G = S — PE — aF (no 
constraint problem). The condition of grand canonical 
stability corresponds to inequality ([M|) for all variations 
bio. By decomposing the perturbation of the eigenmodes 
of the Laplacian, it is easy to show that the system is 
a maximum of grand potential iff /3 > /3ii (where is 
the largest eigenvalue of the Laplacian). This is closely 
related to the Arnol'd theorem (indeed, the grand po- 
tential is equivalent to the Arnol'd energy-Casimir func- 
tional [35| : furthermore, for a linear uj — ip relationship, 
the Arnol'd theorem, which usually provides only a suffi- 
cient condition of grand canonical stability, now provides 
a necessary and sufficient condition of grand canonical 
stability). Since grand canonical stability implies micro- 
canonical stability (but not the converse) we conclude 
that, if /? > /3ii, the system is a maximum of entropy at 
fixed circulation and energy. 

Considering the canonical ensemble, we have to max- 
imize the free energy J = S — f3E at fixed circulation 
(one constraint problem). The condition of canonical 
stability corresponds to inequality ([M)) for all variations 
Slj that conserve circulation. By carefully taking into 
account the constraint on the circulation, Venaille & 
Bouchet (ZTj show that the system is a maximum of free 
energy iff /3 > max{/32i, /?*}. In particular, the states 
with max{/32i, /?*} < P < Pii are stable in the canonical 
ensemble but unstable in the grand canonical ensemble. 
Thus, canonical and grand canonical ensembles are in- 
equivalent [4lj. On the other hand, since canonical sta- 
bility implies microcanonical stability (but not the con- 
verse) we conclude that, if /? > max{/32i, /?*}, the system 
is a maximum of entropy at fixed circulation and energy. 
In particular, the states with E > T'^/{2Ali) = i?ii(r) 
are stable in the canonical ensemble but unstable in the 
grand canonical ensemble [4l| . Note that the states with 
/3 < /3* are unstable in the canonical ensemble (they are 
saddle points of free energy at fixed circulation). This 
result can be obtained directly by considering a pertur- 
bation of the form 6uj = 1 — like in (ii). For this 
perturbation, ST = 0. On the other hand, in the canon- 
ical ensemble, we do not need to impose SE = so that 



this perturbation can be applied to any state leading to 

S'J^l{l3*-m<l>*)~P*{<Pl))>0, (96) 

which proves the result. By contrast, this argument does 
not work in the microcanonical ensemble since the cho- 
sen perturbation does not satisfy SE = for all states. 
Therefore, when /3* > /32i, the states with P21 < /3 < /3* 
are unstable in the canonical ensemble (they are saddle 
points of free energy at fixed circulation) while they are 
metastable in the microcanonical ensemble (they are lo- 
cal maxima of entropy at fixed circulation and energy). 
This is an interesting notion of ensemble inequivalence 
which affects metastable states (Venaille & Bouchet [4l| 
show that the microcanonical and canonical ensembles 
are equivalent for the fully stable states but the case 
of metastable states is not considered in their study). 
In particular, we note that the metastable states with 
/321 </?</?* have negative specific heats (see Fig. [Ij. 
This is allowed in the microcanonical ensemble but not 
in the canonical ensemble. Interestingly, this is the first 
observation of negative specific heats in that context. 

Combining all these results, we conclude that in the 
microcanonical ensemble: 

(a) If /321 < the states are stable for /?>/?*, 
unstable for /3 < /32i and metastable for /32i < /? < /3*, as 
shown in Fig. [T] Therefore, the direct monopole is stable 
for any A^ and the reversed monopole is metastable for 
A2 < (A^O^ 

(b) If P21 > P*- the states are stable for /3 > P21 and 
unstable for /? < /32i, as shown in Fig. [2] Therefore, the 
dipole is stable for A^ < (Ajj^)^ and the direct monopole 
is stable for A^ > {A2i)^. There is no metastable state 
in that case. 



E. The chemical potential 



In Sec. IV C[ we have represented the inverse temper- 
ature /? as a function of A. We shall now study how the 
chemical potential a depends on A. The chemical poten- 
tial is given by a = — F. For the eigenmodes. 



0. 



(97) 



For the solutions of the continuum, assuming T ^ 0, and 
using Eq. dM]), we get a = r/(/3(0) - 1) = T/F{l3). 
Therefore, 



A 



V2E F{P) ■ 

For r = 0, using Eq. ([57]) . we obtain 

a 1 
^= = ± 



(98) 



(99) 



which is a limit case of Eq. ([M)) . The normalized chem- 
ical potential a/y/2E is plotted as a function of A in 
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Figs. |3] and m for a square domain and for a rectangular 
domain of aspect ratio 2, respectively. To plot this curve, 
we have used Eqs. (f85|) and (|98|) . For a given value of 
/?, we can determine A by Eq. ((85|) and a/y/2E by Eq. 

Therefore, we can obtain a/^/2E as a function of 
A parameterized by j3 for the solutions of the continuum. 
For the mixed solutions, P = Pmn is fixed and a/\/2E is 
a linear function of A given by Eq. (|98p . 
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libria is parameterized by /3. The branches corresponding 
to /? > are stable in the grand canonical, canonical 
and microcanonical ensembles and the branches corre- 
sponding to /3 > max{/32i, /?*} are stable in the canonical 
and microcanonical ensembles. The part of the branches 
corresponding to max{/32i, /?»} < /? < arc stable in 
the canonical and microcanonical ensembles but not in 
the grand canonical ensemble. For r < Tc {P21 < /?*), 
the part of the branches corresponding to /32i < /3 < /3* 
are metastable in the microcanonical ensemble and un- 
stable in the other ensembles. 

Remark: in the grand canonical ensemble, the control 
parameter is the chemical potential a and the conjugated 
variable is the circulation F. We must therefore consider 
T{a) by rotating the curves of Figs. [3] and |4] by 90°. 
Only the part of the curve with /? > /3ii (NW and SE 
quadrants) are stable in the grand canonical ensemble. 
There is a first order grand canonical phase transition at 
a = marked by the discontinuity of the circulation T{a) 
between F = ± A"]^ \/2E. Note that there is no metastable 
states in the grand canonical ensemble because the states 
with /3 < Pii are all unstable. 



FIG. 3: Relationship between a/V2E and A in a square do- 
main (case P21 < The solid lines correspond to the sta- 
ble {P > j3*) and metastable (/32i < /3 < /?*) states. Un- 
stable states (/3 < /32i) are represented by the dashed lines. 
The straight lines represent the mixed solutions with constant 
temperature: /3 = /3i2 = /32i, /3 = /322, /? = /323- 
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FIG. 4: Relationship between a/y2E and A in a rectangu- 
lar domain of aspect ratio 2 (case /32i > The solid lines 
correspond to the stable states (/3 > /32i). Unstable states 
(/? < /32i) are represented by the dashed lines. The straight 
lines represent the mixed solutions with constant tempera- 
ture: P = P2I, P ^ P12, P = P22. 

In Figs. [3] and 131 we have represented the series of equi- 
libria containing all the critical points of entropy. If we 
continue the series of equilibria to more and more un- 
stable states, the curve rolls up several times around the 
origin (not shown). As indicated above, the series of cqui- 



F. Description of phase transitions 

We briefiy discuss the nature of phase transitions asso- 
ciated with the maximization problern ([701) and confirm 
the results that Venaille & Bouchet 141[ obtained by a 
different method. We also give a special attention to the 
metastable states that are not considered in (4l| . 

We shall describe successively the caloric curve P{E) 
for a fixed F and the chemical potential curve q;(F) for a 
fixed E. As first observed by Chavanis & Sommeria [26| . 
the nature of the phase transitions depends on the value 
of maxj/?*, /32i}- In a rectangular domain, this quantity 
is determined by the value of the aspect ratio r. We 
must therefore consider two cases successively: r < Tc 
and r > Tc- 



1. Caloric curve 

The caloric curve corresponds to the stable part of the 
series of equilibria I3{E) containing global (stable) and 
local (metastable) maximum entropy states at fixed E 
and F. 

• Let us first consider r < Tc corresponding to 
max{^,,^2i} = P* as in Fig. [H For F = 0, the 
maximum entropy state is the monopole and the caloric 
curve is simply a straight line /3(i?, F = 0) = For 
each value of the energy, we have two solutions with the 
same inverse temperature /3* but different values of the 
chemical potential a(F = 0,i?) = ±ao (see Fig. 
One solution is a monopole with positive vorticity at 
the center (MP) and the other solution is a monopole 
with negative vorticity at the center (MN). For F = 0, 
these solutions have the same entropy. Thus, the branch 
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l3{E,r = 0) = /3* is degenerate. For T ^ 0, the 
caloric curve /?(£', F 7^ 0) can be deduced easily from 
Fig. [T|63I |. The global maximum entropy state is the 
direct monopole (for F > the vorticity is positive at 
the center (MP); for F < the vorticity is negative at 
the center (MN)) for any E. For E > £'^i(F) = 2(i^' 

the reversed monopole is metastable (local entropy max- 
imum). Note that the metastable states have negative 
specific heats C = dE/d{l/(3) < 0. The caloric curve 
/3{E) does not present any phase transition. 

• Let us now consider r > Tc corresponding to 
max{/3*,;32i} = /?2i as in Fig. El For F = 0, the 
maximum entropy state is the dipole and the caloric 
curve is simply a straight line /3{E,r = 0) = /32i. For 
each value of the energy, we have two solutions with the 
same inverse temperature and the same chemical poten- 
tial a(F = 0, i?) = (see Fig. [4]). One solution is a 
dipole (-I-, — ) with positive vorticity on the left and the 
other solution is a dipole (—,-!-) with negative vortic- 
ity on the left (in Fig. [51 we have only represented the 
dipole (— , -I-)). For F = 0, these solutions have the same 
entropy. Thus, the branch f3{E,r = 0) = /?2i is degen- 
erate. For F 7^ 0, the caloric curve I3{E,T ^ 0) can 
be deduced easily from Fig. [51 The maximum entropy 
state is the asymmetric (mixed) dipole (-1-,—) or (—,-!-) 
for E > i?2i(F) = 2(A\y^ ^-"^"^ ^^"^ direct monopole for 
E < i?2i(r) = 2(1'' )3 (^^^ reversed monopoles are un- 
stable). The caloric curve f3{E) presents a second order 
phase transition marked by the discontinuity of ^{E) 
at E = i;^i(F). 



2. Chemical potential curve 

The chemical potential curve corresponds to the sta- 
ble part of the series of equilibria a(F) containing global 
(stable) and local (metastable) maximum entropy states 
at fixed E and F. 

• Let us first consider r < Tc corresponding to 
max{/3*,^2i} = /5* as in Fig. [3l The global maximum 
entropy state is the monopole for any value of F. Con- 
sidering only fully stable states (global entropy maxima), 
there is a first order phase transition at F = marked by 
the discontinuity of Q;(r) while the entropy is continuous. 
When we pass from positive F to negative F, we pass dis- 
continuously (in terms of a but not in terms of /? or S) 
from the monopole (MP) to the monopole (MN). In fact, 
due to the presence of long-lived metastable states (see 
Sec. IVI Cp . we remain in practice on the monopole (MP) 
until the metastable branch disappears. Then we jump 
on the monopole (MN) with discontinuity of a (and /3 
and S). This corresponds to a zeroth order phase tran- 
sition. 

• Let us now consider r > Tc corresponding to 
max{/3*, /32i} = /32i as in Fig. [H The global max- 
imum entropy state is the asymmetric (mixed) dipole 



for |F| < F2i(-E') = V2EA21 and the direct monopole 
for |r| > r2i{E). There are two second order phase 
transitions marked by the discontinuity of ^ (F) at F = 

±r2,{E). 

3. Phase diagram 

The phase diagram in the (t. A) plane, including the 
metastable states, is plotted in Fig. [5l Depending on 
the values of A and r (and depending on the history 
of the system in the zone of metastability) , the maxi- 
mum entropy state is a dipole (D), a monopole (MP) or 
a monopole (MN). If we fix the circulation F, we obtain 
the phase diagram in the {t,E) plane. For F 7^ 0, it 
shows the appearance of a second order phase transition 
in /?(£') for t > Tc (for F = there is no phase transition). 
If we fix the energy E, we obtain the phase diagram in 
the (r, F) plane. As noted by Venaille & Bouchet [4]| . 
the point (F = 0, r = t^) is a bicritical point marking 
the change from a first order to two second order phase 
transitions in a(F). 
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FIG. 5: Phase diagram in the (r. A) plane showing the domain 
of stability of the direct monopoles and dipoles. We have 
indicated by a dashed line the domain of metastability of the 
reversed monopoles. 

Remark: for illustration, we have described the phase 
transitions in the case of a rectangular domain and for the 
Laplacian operator. The generalization to an arbitrary 
domain and a linear operator L is straightforward. In 
that case (32i is replaced by f3[ (the first eigenvalue of 
L with zero mean) and /3ii is replaced by /3" (the first 
eigenvalue of L with non zero mean) . 



VI. RELAXATION TOWARDS MINIMUM 
ENSTROPHY STATES 

We shall now illustrate numerically the phase transi- 
tions discussed previously using the relaxation equations 
introduced in Appendix [D1 These relaxation equations 
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can serve as numerical algorithms to compute maximum 
entropy states or minimum enstrophy states with rele- 
vant constraints. Their study is also interesting in its 
own right since these equations constitute non trivial dy- 
namical systems leading to rich bifurcations. Although 
these relaxation equations do not provide a parametriza- 
tion of 2D turbulence (we have no rigorous argument for 
that), they may however give an idea of the true evolution 
of the flow towards equilibrium. In that respect, it would 
be interesting to compare these relaxation equations with 
large eddy simulations (LES) of 2D turbulence. This will, 
however, not be attempted in the present paper. 



A. Relaxation equations 

We shall numerically solve the relaxation equation of 
Sec. ID 31 For simplicity, we shall ignore the advective 
term since we are just interested in describing the bifur- 
cations between the different equilibrium states. Then, 
by a proper rcscaling of time, we can take D — I with- 
out loss of generality. The relaxation equation (|D43P 
becomes 



(100) 



with the boundary condition u = —a{t) on the edge 
of the domain. The Lagrange multipliers /3(t) and a(t) 
evolve in time according to Eqs. (|D44p and (jD45[) in or- 
der to conserve the circulation and the energy. This leads 
to 



m 



T{il;) - 2E 

m - W 



(101) 



temperature /?* . It is the maximum entropy state at fixed 
circulation and energy. In fact, for F = 0, this solution 
is degenerate since the monopoles (MP) and (MN) have 
the same entropy. 

Let us confront these theoretical results to a direct nu- 
merical simulation of Eq. (jlOOp . Starting from a generic 
initial condition (made of Gaussian peaks with positive 
and negative vorticity symmetrically distributed in the 
domain to assure F = 0), we numerically find that the 
system spontaneously relaxes towards the dipole and re- 
mains in that state for a long time (see Fig. [5]) although 
this state is predicted to be unstable (see Sec. IV C[) . This 
simple numerical experiment shows that unstable states 
can be long-lived. In fact, the dipole is a saddle point of 
entropy so that it is unstable only for very specific pertur- 
bations. If these perturbations arc not generated spon- 
taneously during the relaxation process, the system can 
remain frozen in a saddle point for a long time. Another 
reason why the dipole has a long lifetime is due to the fact 
that the entropies of the monopole (stable) and dipole 
(unstable) are very close for F = since —46.5 
and 1321 ~ —49.3. To check that the dipole is really 
unstable, we have introduced by hands (see the arrow in 
Fig.|6l) an optimal perturbation of the form Suj = 1-/3*0* 
(see Sec. IVDp . In that case, the dipole is immediately 
destabilized and the system quickly relaxes towards the 
monopole which is the maximum entropy state in that 
case. In the case shown in Fig. [6l we obtain a monopole 
(MP). If we introduce an optimal perturbation with the 
opposite sign, wc get the monopole (MN). If we do not 
introduce any perturbation by hand and just let the sys- 
tem evolve with the numerical noise, the dipole finally 
destabilizes but this takes a long time (not shown) of the 
order t ^ 400. 



The rate of increase of entropy (neg-enstrophy) is 



S= {uj + (3i' + ay dr > 0. 



(102) 



(103) 



Therefore, the relaxation equation (jlOOp with the con- 
straints (jlOip and (|102p relaxes towards the maximum 
entropy state at fixed circulation and energy. Saddle 
points of entropy are linearly unstable to some pertur- 
bations (in particular those described in Sec. IV D[) . 




B. Geometry induced phase transitions and 
persistence of saddle points 

We first consider the case of a square domain (r = 1 < 
Tc) and take F = 0. For these values of parameters, the 
relaxation equation (|100p admits an infinite number of 
steady states that are the solutions of Eq. ((75)) . How- 
ever, the only stable solution is the monopole with inverse 



FIG. 6: Starting from a generic initial condition with F = 
in a square domain, the system relaxes towards a dipole (first 
plateau) although this solution is unstable (saddle point). At 
t — 20 (see arrow), an optimal perturbation is applied to the 
dipole which quickly destabilizes in a stable monopole (second 
plateau). In the absence of optimal perturbation, the system 
can remain frozen in the dipole for a long time. 

We now consider a rectangular domain with aspect ra- 
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tio r = 2 > Tc and again take F = 0. In that case, the 
maximum entropy state at fixed circulation and energy is 
the dipole and the monopole is unstable (saddle point). 

Starting from a generic initial condition, the system 
spontaneously relaxes towards the dipole and remains in 
this state even if very large perturbations are applied (not 
shown). By contrast, if we start from the monopole, we 
numerically observe that the system remains in that state 
for a very long time although this state is unstable (see 
Sec. IV Cp . If we apply by hands (see the arrow in Fig. [7]) 
an optimal perturbation of the form Suj = ^21 (see Sec. 
IV D[) . the monopole is immediately destabilized and the 
system quickly evolves towards the dipole (Fig. [7]) which 
is the maximum entropy state in that case. In the ab- 
sence of applied perturbation, we have not observed the 
destabilization of the monopole on the timescale achieved 
in the numerical experiment (however, if we add the ad- 
vection term, the dipole is formed on a time of the order 
t ~ 400). 

-35 1 , , , > , , , 1 




.55 1 1 1 1 \ 1 \ 1 

5 10 15 20 25 30 35 40 

t 

FIG. 7: Starting from a monopole with F = in a rectangu- 
lar domain (r > Tc), the system remains in that state for a 
long time (first plateau) although this state is unstable (sad- 
dle point). At t = 10 (see arrow), an optimal perturbation 
is applied to the monopole which quickly relaxes towards a 
dipole (second plateau). In the absence of optimal perturba- 
tion, the system can remain frozen in the monopole state for 
a long time. 

In conclusion, this numerical study reveals that even 
unstable states (saddle points of entropy) can be natu- 
rally selected by the system and persist for a long time. 
Indeed, these states are destabilized by a very partic- 
ular type of perturbations (that we call optimal) and 
such perturbations may not be necessarily generated by 
the internal dynamics of the system. This suggests that 
the system can be frozen for a long time in a quasi sta- 
tionary state (QSS) that is not necessarily a stable or 
metastable steady state of the 2D Euler equation. It 
can even be an unstable saddle point! This observation 
has been made on the basis of the relaxation equations 
that are constructed so as to relax towards a maximum 
entropy state. However, the same phenomenon could 
appear for real flows described by the Eulcr or Navier- 



Stokes equations in numerical simulations and laboratory 
experiments. This could be interesting to study in more 
detail. 



C. Metastability and hysteresis 

We shall now describe the hysteretic cycle predicted 
by statistical mechanics (based on the neg-enstrophy) in 
a domain with aspect ratio r < Tc. In Fig. [51 we plot 
the entropy S/ E as a function of the control parameter 
A. We shall assume that the energy E is fixed so that 
A basically represents the circulation F. The hysteresis 
is due to the presence of metastable states (local entropy 
maxima) when -A21 < A < A21. For < A < A21, the 
global maximum entropy state is the direct monopole 
(MP) while the reversed monopole (MN) is metastable. 
For —A21 < A < 0, the global maximum entropy state is 
the direct monopole (MN) while the reversed monopole 
(MP) is metastable. Depending on how it has been pre- 
pared initially, the system can be found in the stable or 
metastable state. 




A 

FIG. 8: S/E ratio as a function of A in a square domain. 

We start from a state with large A corresponding to 
positive temperature (/3 > 0). In that case, the positive 
vorticity has the tendency to accumulate on the bound- 
ary of the domain. If we reduce A, we enter in the region 
of negative temperature states {(3 < 0). In that case, 
the positive vorticity has the tendency to accumulate at 
the center of the domain. For A > 0, the (global) maxi- 
mum entropy state is the monopole (MP). For A = 0, we 
expect a first order phase transition from the monopole 
(MP) to the monopole (MN) (see Sec. IVC)) marked by 
the discontinuity of the chemical potential a (while /3 and 
S are still continuous). In fact, for —A21 < A < 0, the 
monopole (MP) is metastable and robust so that the sys- 
tem remains on this branch. Therefore, in practice, the 
first order phase transition does not take place. How- 
ever, for A < —A'21, the branch of monopoles (MP) be- 
comes unstable and the system jumps to the branch of 
direct monopoles (MN) which correspond to global en- 
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tropy maxima. This is marked by a discontinuity of en- 
tropy (zeroth order phase transition). If we decrease A 
sufficiently, we enter in the region of positive tempera- 
ture states (/? > 0). In that case, the negative vorticity 
has the tendency to accumulate on the boundary of the 
domain. If we now increase A the system follows the 
branch of monopoles (MN) which is stable for A < and 
metastable for < A < Aj^. Again, the first order phase 
transition at F = does not take place. For A > A21, 
the branch of monopoles (MN) becomes unstable and the 
system jumps to the branch of direct monopoles (MP) 
which correspond to global entropy maxima. We have 
thus followed an hysteretic cycle as illustrated in Figs. [8] 
and[9l 




FIG. 9: Hysteretic cycle in a square domain, obtained by nu- 
merical integration of Eq. (|100|l . We have represented A 
(black) and /3 (red) as a function of time. Starting from 
a stable state with A G [0; A21] (MP), the system is regu- 
larly perturbed: at t = 10,250,490,730,970,1210, we add 
to the vorticity distribution the sum of a negative Gaus- 
sian peak and an eigenmode 0.1'i/'2i, and let the system re- 
lax. The effect of the Gaussian peak is to decrease A, 
while the eigenmode destabilizes the unstable states. For 
< A < A21 we follow the stable branch (MP) of Fig. [S] and 
for — A21 < A < 0, we follow the metastable branch (MP). 
For A < — A21, the metastable solutions (MP) no longer ex- 
ist, and the system jumps to the upper branch (MN) of Fig. 
m At f = 1450, 1690, 1930, 2170, 2410, we add to the vortic- 
ity distribution the sum of a positive Gaussian peak and an 
eigenmode 0.05V'2i ■ The value of A is then increased, and 
we follow the stable branch (MN) for A G [— A2i;0] and the 
metastable branch (MN) for A G [0; A21]. When A > A21, the 
metastable solutions (MN) no longer exists, and the system 
jumps to the upper branch (MP) of Fig. [S] 



D. Bifurcations in the presence of a noise 

For F = in a square domain, the monopoles (MP) 
and (MN) are stable and have the same entropy but re- 
main quite distinct states (with opposite velocity). This 
corresponds to a parity breaking for the final organiza- 
tion of the system |26j . In the presence of forcing, we 



expect to observe random transitions between these two 
solutions |64[ similar to those observed experimentally by 
Sommeria [43| for 2D turbulence forced at small scale in 
a square box. Indeed, we are in a situation similar to the 
case of a bistable system. To observe such transitions, 
one possibility is to introduce a stochastic noise in the 
relaxation equation (jlOOp . Unfortunately, for a simple 
white noise, we did not observe any transition and we 
have not been able to find the properties of forcing that 
allow such transitions to appear. This may be due to the 
high entropic barriere created by the unstable (dipole) 
solution. Therefore, in order to illustrate the main idea, 
we shall introduce a simple effective model. 

The relevant order parameter is the chemical potential 
a which takes the values iao for the (stable) monopoles 
(MP) and (MN) and the value a = for the (unsta- 
ble) dipole (see Fig. [S]). We shall now introduce an 
entropic function S{a) modeled by a symmetric func- 
tion with three bumps (two maxima and one minimum). 
Since we know the entropy (by unit of energy) of the 
monopoles Smonopoies = P* and the entropy of the dipole 
Sdipoie = hi , we find that 



S{a) = (/321 - /?*) 



1-1^ 



n 2 



(104) 



When a forcing is present, we can propose that a becomes 
a stochastic variable described by a Langevin equation of 
the form 



^=M5'(a) + 



(105) 



where •q{t) is a white noise. In the absence of forcing, 
Eq. (|105|) relaxes towards one maximum of S{a), the 
monopole (MP) or the monopole (MN), and stay there 
permanently. In the presence of forcing, Eq. (|105p de- 
scribes random transitions between these two states (see 
Fig. [TU)) . This is the classical bistable system that has 
been studied at length in statistical mechanics and Brow- 
nian theory [soj . 

Random transitions have been observed in various 
physical systems in fluid mechanics (see, e.g., ^6, 49,, Sj 
and references therein). In the present study, we have 
considered random transitions between a monopole (MP) 
and a reversed monopole (MN). They are associated with 
the first order phase transition that takes place in a 
square domain when the QSS has a linear uj — ip rela- 
tionship. It would be interesting to see if they can be ob- 
tained directly from the forced Navier-Stokes equations 
in situations where the co — ip relationship is close to lin- 
ear. Random transitions between a unidirectional flow 
and a dipole have been obtained recently by Bouchet & 
Simonnet [46t by solving numerically the forced Navier- 
Stokes equations in periodic domain. However, the situ- 
ation is different (and more complex) because these two 
states are characterized by different uj — ip relationships. 
Indeed, the forcing can change the shape of oj{ip). In the 
situation that we consider, the shape oiojiip) remains the 
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FIG. 10: Solution of the stochastic equation (|105p for fi = 
1.0 and D — 1.25 showing random transitions between the 
monopoles (MP) and (MN). The dipole is always unstable. 



same (linear) but the equation Aip ~ —uj{ip) determining 
the QSS can admit two stable solutions (MP) and (MN). 
This situations is closer to that of a bistable system and 
would be interesting to study numerically. 



VII. CONCLUSION 

In this paper, we have studied the maximization of 
the Miller-Robert-Sommeria entropy Smrs at fixed en- 
ergy E, circulation F and microscopic cnstrophy T2^' 
and proved the equivalence with the minimization of the 
macroscopic enstrophy Fj ^' at fixed energy E and cir- 
culation F. This provides a justification of the minimum 
enstrophy principle from statistical theory when only the 
microscopic enstrophy is conserved among all the Casimir 
invariants. We have suggested that relevant constraints 
(such as the microscopic enstrophy) are selected by the 
properties of forcing and dissipation. Our simplified ther- 
modynamic approach leads to a mean flow characterized 
by a linear uj — ip relationship and Gaussian fluctuations 
around it. Such states can be relevant to describe certain 
oceanic flows [6l-l20|. More general flows with nonlinear 
U — tjj relationships (and more general fluctuations) can 
be constructed in principle by keeping other Casimir con- 
straints in addition to the microscopic enstrophy. 

We have studied the minimization of enstrophy at fixed 
energy and circulation and analyzed the corresponding 
phase transitions with the approach of Chavanis & Som- 
meria [1^ . We have discussed the link with the approach 
of Venaille & Bouchet [4l| . We have proposed relaxation 
equations to solve this minimization problem (see [35| 
for generalizations) and used them to illustrate the phase 
transitions. 

One interesting result of the simulations is the observa- 
tion that saddle points of entropy can be relevant in the 
dynamics. Indeed, these states are unstable only for par- 



ticular perturbations that are not necessarily generated 
spontaneously by the system. As a result, they can be 
long-lived and robust. This observation may have inter- 
esting application in the case of von Karman flows since 
it is found that Beltrami states are saddlepoints of en- 
ergy at fixed helicity, not energy minima [38j. Still, it is 
observed experimentally [s^ [s^ that they are long-lived 
and robust. 

We have also discussed in detail the metastable states 
that were not considered in the study of Venaille & 
Bouchet [4l|. For long-range interactions, metastable 
states (local entropy maxima) are long-lived and they 
are as much important as fully stable states (global en- 
tropy maxima). Interestingly, these metastable states 
have negative specific heats leading to a form of ensem- 
ble inequivalence between microcanonical and canonical 
ensembles (while these ensembles are equivalent at the 
level of fully stable states [4l|). These metastable states 
can lead to an hysteresis and to random transitions be- 
tween direct monopoles and reversed monopoles. Such 
transitions can also arise in more realistic fluid systems 
and can have some importance in oceanography and me- 
teorology giliilli- 

A last remark may be in order. The MRS statistical 
theory of the 2D Euler equation, which is the most ba- 
sic and the most rigorous, takes into account an infinite 
number of constraints. When applied to real flows, this 
is clearly unphysical and this leads to practical difficul- 
ties. It has been a subject of intense debate in the last 20 
years to find a practical way to deal with the constraints. 
Different approaches have been proposed: some consider 
a point vortex approximation where only the energy and 
the number of vortices in each species matter [2^ . some 
consider since the start only a finite number of inviscid 
constraints [l8l - l20| . some consider a strong mixing (or 
low energy) limit of the MRS statistical theory which 
makes a hierarchy among the Casimir constraints [26j . 
and some model the vorticity fluctuations by a prior dis- 
tribution [29l - l3T| . In our recent works [H, [s^] , including 
the present one, we have not tried to determine which 
approach, if any, is the "best" . For the moment, we 
just present different ways to deal with the constraints 
and systematically study the corresponding variational 
principles. We have also extended these variational prin- 
ciples to 3D axisymmctric flows [s^. These variational 
principles have a long history in 2D turbulence and MHD 
and one virtue of our papers is to put several variational 
principles in correspondance. The determination of the 
"best" approach is still a matter of debate and research. 



Appendix A: Equivalence between (|29p and (|54p 

In Sec. llVi we have shown the equivalence of and 
(j46p for global maximization. In this Appendix, we show 
the equivalence of and pS)) for local maximization, 
i.e. p{r,a) is a (local) maximum of S[p\ at fixed E, F, 



and normalization if, and only if, the correspond- 
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ing coarse-grained vorticity CL'(r) is a (local) minimum 
of r2^ [aJ] at fixed E and T. To that purpose, we show 



the equivalence between the stability criteria (f29| and 
(|5^. We use a general method similar to the one used 
in j32l. [521 . [53! in related problems. 

We shall determine the optimal perturbation <5p*(r, cr) 
that maximizes S'^J[Sp] given by Eq. (f29| with the con- 
straints (5ZJ J (5/9cr da, dT2^' = J Spa^ dadv = and 
J Spda ~ 0, where StJ(r) is prescribed (it is only as- 
cribed to conserve circulation and energy at first order). 
Since the specification of 5uj determines Sip, hence the 
second integral in Eq. (|29p . we can write the variational 
problem in the form 



S —. 



-pS (^J Spa^ dcrdr^ - J C(r)(5 (^J Spda^ dr = 0, (Al) 

where A(r), p and C(r) are Lagrange multipliers. This 
gives 

Sp4r, a) = -p(r, a){pa^ + X{v)a + C(r)), (A2) 

and it is a global maximum of (5^J[(5p] with the previ- 
ous constraints since S^{6'^J) = — J ^^^^^^ drda < 
(the constraints are linear in Sp so their second varia- 
tions vanish). The Lagrange multipliers are determined 
from the above-mentioned constraints. The constraints 
J Spda = and SuJ = J Spa da lead to 



C(r) + A(r)aJ(r)-F/^w2(r) =0, 



(A3) 



C(r)w(r) + A(r)w2(r) + puj^{r) = -SuJ{r). (A4) 

Now, the state p(r, a) corresponds to the Gaussian distri- 
bution ([77]) . Therefore, we have the well-known relations 
uj^{r) = a7^(r) + uj2 and uj-^{r) = ZU^(r) -I- 3UJ{r)u!2 where 
LU2 = is uniform. Substituting these relations in Eqs. 
[3]) and (jA4[) . and solving for A(r) and C(r), we obtain 



A(r 



Soj{r) 

U)2 



2/iw(r), 



C(r) = ^^(5aJ(r) + puJ^ir) - puj2- 



1^2 



(A5) 



(A6) 



Therefore, the optimal perturbation (|A2I) can be rewrit- 
ten 



Sp* 



Suj 



ijJ2 



((T - Cli) + ^ { ((T - Cj)^ - W2 } 



(A7) 



The Lagrange multiplier p is determined by substituting 
this expression in the constraint J Spa"^ drda = 0. Using 
the well-known identity u}'^{r) = aJ*(r) -I- 6a;2a7^(r) + 3aj| 



valid for a Gaussian distribution, we obtain after some 
simplifications 



M = 



J loSlo dr 



(A8) 



Therefore, the optimal perturbation (jA2p is given by Eq. 
(|Tf)) with Eq. Since this perturbations maximizes 

S'^J[Sp] with the above-mentioned constraints, we have 
S^J[Sp] < S^J[Sp*]. Explicating (5^ J using Eqs. (|T7)) 
and (jA8[) . we obtain after simple calculations 



{Sujf dr 



jSlu dr 



^-/3 I SuSipdr. (A9) 



The r.h.s. returns the functional appearing in Eq. (|43p . 
We have already explained in Sec. IIVBI that for the 
class of perturbations that we consider [ST = SE = 0) 
the second integral vanishes. Therefore, the foregoing 
inequality can be rewritten 



S'^ASp] < - 



2UJ2 



{5Zjfdr--l3 I SZJSiPdr, (AlO) 



where the r.h.s. is precisely the functional appearing in 
Eq. ([M]). Furthermore, there is equality in Eq. (|A10|) iff 
Sp = Sp^,. This proves that the stability criteria and 
([54]) are equivalent. Indeed: (i) if inequality ([54)) is ful- 
filled for all perturbations Suj that conserves circulation 
and energy at first order, then according to Eq. (jAlO[) . 
we know that inequality is fulfilled for all pertur- 
bation Sp that conserves circulation, energy, fine-grained 
enstrophy and normalization at first order; (ii) if there 
exists a perturbation SZJ* that makes (5^ J[5ZU] > 0, then 
the perturbation Sp* given by Eq. (|A7|) with Eq. (|A8|) 
and SlJ = SZJ* makes S'^J[Sp] > 0. In conclusion, the 
stability criteria (|29| and (|54|) are equivalent. 



Appendix B: Eigenvalues and eigenfunctions of the 
Laplacian in a rectangular domain 

We define the eigenfunctions and eigenvalues of the 
Laplacian by 



(Bl) 



with -iAn = on the domain boundary. These eigenfunc- 
tions are orthogonal and normalized so that (V'nV'm) = 
Snm- Since — /(V?/'„)^dr = /3„ J ^^dr, we note that 
Pn < 0. Following Chavanis & Sommeria [2^, we dis- 
tinguish two types of eigenmodes: the odd eigenmodes 
-0^ such that {ip'^) = and the even eigenmodes V'" such 
that (-0") 7^ 0. We note /3'j and /3" the corresponding 
eigenvalues. 

In a rectangular domain of unit area whose sides are 
denoted a — -y/r and b = 1/\/t (where t = a/b is the 
aspect ratio), the eigenmodes and eigenvalues are 



ipmn = 2 sm{rmrx / \/t) sin(n7r-\/r2/). 



(B2) 
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(B3) 



where the origin of the Cartesian frame is taken at the 
lower left corner of the domain. The integer m > 1 gives 
the number of vortices along the a;-axis and n > 1 the 
number of vortices along the y-axis. We have (ipmn) = 
if m or n is even and {tpmn) ^ if m and n are odd. 

The differential equation (|82|) can be solved analyt- 
ically by decomposing the field on the eigenmodes 



as (j) = J2 
obtain 



'4>mn and using the identity 1 



This yields Eq. (|83p from which we 



We note in particular that 



dl3 ' 



I M\ a I Mix Pmn{'4'mn) „ 



(B4) 



(B5) 



(B6) 



Appendix C: Temporal evolution of the different 
modes 



The relaxation equation (jlOOp can be solved analyti- 
cally by decomposing the vorticity and the stream func- 
tion on the eigenmodes of the Laplacian. Using the 
Poisson equation, we get oj{r,t) = ""(*)V'ri(r) and 
= E„^n(i)^«(i-) with hn{t) = -an{t)IPn- Sub- 
stituting these expressions in Eq. (|100p and using the 
identity 1 = X]n(V'n)'0n, we obtain the ordinary differen- 
tial equations 

^+(l-^)a« = -«W(V'.), (CI) 

for all n. The evolution of the Lagrange multipliers is 
given by Eqs. (UnB and (HUH) with {iP) = J2nbn{t){i^n) 
and (tp^) = X]n^n(*)- The modes are coupled through 
the Lagrange multipliers in order to assure the conserva- 
tion of energy and circulation. 

In the grand canonical description in which /3 and a 
are constants, the foregoing differential equation can be 
integrated straightforwardly, yielding 



a„(0 = a„(0)-F 



- - P/Pn 



-(l-/3//3„)t 



1 - P/Pn ■ 



(C2) 



In that case, a steady state of the relaxation equation 
is stable iS f3 > /3[ where P[ is the largest eigenvalue of 
the Laplacian. The condition /3 > is a necessary and 
sufficient condition for the steady state to be a global 



maximum of the grand potential G = S — PE — aV. 
That functional is related to the Arnol'd energy- Casimir 
functional used to settle the nonlinear dynamical stability 
of a steady state of the 2D Euler equation [3^1 ■ 



Appendix D: Relaxation equations 

1. Relaxation equations associated with the 
maximization problem (|16p 

In this Appendix, we construct relaxation equations as- 
sociated with the maximization problem (|16|) correspond- 
ing to the energy-enstrophy-circulation statistical theory. 
These relaxation equations can serve as a numerical al- 
gorithm to solve this constrained maximization problem. 
In the past, Robert & Sommeria [s^l have proposed re- 
laxation equations that conserve all the Casimirs and in- 
crease the entropy. Here, we use a different approach 
because we want to conserve only the microscopic en- 
strophy (not all the Casimirs). Thus, the form of the 
relaxation equations will be different. In particular, they 
will involve a current in the space of vorticity levels a 
[3^ . [sst instead of a current in the space of positions [sj] . 

We construct a set of relaxation equations that increase 
S[p\ while conserving i?, F and T2^' using a Maximum 
Entropy Production Principle. The dynamical equation 
that we consider can be written as 



dp 
dt 



u-Vp 



dJ_ 



(Dl) 



where J is an unknown current to be chosen so as to 
increase S[p\ while conserving the constraints. The local 
normalization J pda = 1 is satisfied provided that J — > 
as cr ±00. Multiplying Eq. (jDl[) by a and integrating 
over the levels, we get 



duj „ 
— + u • VZU 

ot 



Jda = X. 



(D2) 



Next, multiplying Eq. (jPip by and integrating over 
the levels, we obtain 



— u • Vw2 = 2 / Jada. (D3) 

From Eqs. ((D2)) and (|D3l) . we find that 

^ -t- u • Vw2 = 2 / J(cr - uj)d(T. (D4) 
at J 

Using Ec^. (jPlI) . the time variations of S[p\ are given by 

J dp 



S 



p da 



drda, 



and the time variations of E, F, T^'^' are given by 



E = / Ji^drda = 0, 



(D5) 



(D6) 
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T = J drda = 0, 



J a drda ~ 0. 



(D7) 



(D8) 



Following the Maximum Entropy Production Principle, 



we maximize S with E = T 
tional constraint 



-pfg- 



2p 



da < C{r,t), 



and the addi- 



(D9) 



putting some physical bound on the diffusion current. 
The variational principle can be written in the form 



5S - l3{t)6E - a{t)St - a2it)Stl-^- 



( [ -T^da 1 dr 



0, 



(DIO) 



where (3(t), a(t), a2{t) and D[r,t) are time dependent 
Lagrange multipliers associated with the constraints. 
This leads to the following optimal current 



J = -D 



dp 



- p{^(t)'iP + a{t) +2a2{t)a) 



(Dll) 



Therefore, the relaxation equation for the vorticity dis- 
tribution is 



dp 



^+p(/3(i)?A + a(i)+2a2(i)tT) 
aa 



Integrating Eq. (|Dlip over cr, we obtain 

X^-B (/3(t)V + oi(t) + 2a2(t)cJ) • 



|(D12) 



(D13) 



Inserting Eq. (|D13p into Eq. ((D2|) leads to the following 
relaxation equation for the mean flow 



'dt 



u • Voj 



-D {I3{t)i; + a{t) + 2a2{t)uJ) . (D14) 



For the boundary condition, we shall take I3{t)ip + a{t) + 
2a2 {t)LJ = on the domain boundary so as to be consis- 
tent with the equilibrium state where this quantity van- 
ishes in the whole domain. Since ip = on the boundary, 
we finally get uj = —a{t)/{2a2it)) on the domain bound- 
ary. A relaxation equation can also be written for the 
centered variance 0J2- Using Eqs. (|D11|) and (|D4[) . wc 
obtain 



^ + u • Vw2 = 2L' (1 - 2a2{t)L02) ■ 
ot 



(D15) 



Finally, in Eqs. (|D12|) . (|D14|) and (|D15|) . the Lagrange 
multipliers evolve so as to satisfy the constraints. Sub- 
stituting Eq. (|DTT|) in Eqs. dHSl), (lD7| and dDS]), we 
obtain the algebraic equations 



(^2)/3(t) + (^)a(t)+4i;a2(i)=0, 



(D16) 



{iP)P{t) + a{t) + 2Ta2{t) = 0, 



2EI5{t) + Tait) + 2T{-^a2{t) = I. 



(D17) 
(D18) 



where (X) = J X dr. Substituting 9p/(9(T taken from Eq. 
(jPlip in Eq. ((D5|) and using the constraints (|D6|) - (|D8|) . 
we easily obtain 



J2 

— — drda, 
Dp 



(D19) 



so that S > provided that D is positive. On the other 
hand 5* = iff J = leading to the Gibbs state (^5]) . 
From Lyapunov's direct method, we conclude that these 
relaxation equations tend to a maximum of entropy at 
fixed energy, circulation and microscopic enstrophy. Note 
that during the relaxation process, the distribution of 
vorticity is not Gaussian but changes with time according 
to Eq. (|D12p . The vorticity distribution is Gaussian only 
at equilibrium. Therefore, these relaxation equations de- 
scribe not only the evolution of the mean flow according 
to Eq. (jD14[) but also the evolution of the full vorticity 
distribution according to Eq. (jD12[) . We stress, however, 
that these equations are purely phenomenological and 
that there is no compelling reason why they should give 
an accurate description of the real dynamics. However, 
they can be used at least as a numerical algorithm to 
compute the statistical equilibrium state. Indeed, these 
equations can only relax towards an entropy maximum at 
fixed energy, circulation and microscopic enstrophy, not 
towards a minimum or a saddle point that are linearly 
unstable with respect to these equations [65 1. 



2. Relaxation equations associated with the 
maximization problem (|36p 

We shall now introduce a set of relaxation equations as- 
sociated with the maximization problem p6p. We write 
the dynamical equation as 

^ u • = X, (D20) 

where X is an unknown quantity to be chosen so as to 
increase S[uj] while conserving E, F and Fj'^'. The time 
variations of S arc given by 

S^-^J:,Xdr, (D21) 

where fl2{t) is determined by the constraint on micro- 
scopic enstrophy leading to 



^2it) = T^-^- - J uJ^dr, 



(D22) 



at each time. On the other hand, the time variations of 
E and F are 



(D23) 
(D24) 



E = J X2jjdr = 0, 
f = Xdr^O. 
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Following the Maximum Entropy Production Princi- 
ple, we maximize S with E = T = (the conservation 
of microscopic enstrophy has been taken into account in 
Eq. (1022^ ) and the additional constraint 



<C(r,t). 



(D25) 



The variational principle can be written in the form 

6S - - ait)6t - J -^^S (^^^ dr - 0, 

(D26) 



and it leads to the optimal quantity 
Inserting Eq. ((D27)) in Eq. ((D20)) . we obtain 



(D27) 



^ + u • VoJ = -Z? (^/3(0 V + a{t) + , (D28) 

with aJ = —a(t)^2(t) on the domain boundary. The La- 
grange multipliers evolve so as to satisfy the constraints. 
Substitu ting E q. (|D27| in Eqs. ((D23| - ((D24)) and recall- 
ing Eq. (|D22p . we obtain the algebraic equations 



- rf-9- _ I ^2 



2E 



(D29) 



(D30) 



(D31) 



Substituting uJ taken from Eq. (|D27|) in Eq. (|D2ip and 
using the constraints (jD23P - (|D24p . we easily obtain 



S 



D 



dr, 



(D32) 



so that iS* > provided that D is positive. On the other 
hand S' = iff X = leading to the condition of equilib- 
rium (j42p . From Lyapunov's direct method, we conclude 
that these relaxation equations tend to a maximum of 
entropy at fixed energy, circulation and microscopic en- 
strophy. 

The relaxation equation (|D28p is similar to Eq. (jD14p 
but the constraints determining the evolution of the La- 
grange multipliers are different. More precisely, Eqs. 
(|D30)) and (|MT1) are equivalent to Eqs. 0T7\ and ((DTH)) 
but Eq. (jPiep has been replaced by Eq. (|D29p . In- 
deed, in the present approach, the vorticity distribution 
is always Gaussian during the dynamical evolution. It is 
given by Eq. p2p at any time, i.e. 



p(r,cr,t) = 



1 



■-TJ{r,t))^ 
2Q2(t) 



(D33) 



By contrast, in the approach of Sec. ID 11 the vorticity 
distribution changes with time. Therefore, the dynamical 
evolution is different. However, in the two approaches, 
the equilibrium state is the same, i.e. it solves the maxi- 
mization problem (|16p. This is sufficient if we use these 
relaxation equations as numerical algorithms to compute 
the maximum entropy state. 

Remark: Using Eqs. (|D20p - (|D2ip . it is easy to show 
that f2^- = -2VL2{t)S so that fa^' < since ^2(1) > 
(by Schwartz inequality). Therefore, the macroscopic en- 
strophy decreases monotonically through the relaxation 
equations. This is to be expected since the maximiza- 
tion problem (j36]) is equivalent to the minimization of 
the macroscopic enstrophy at fixed energy and circula- 
tion (see Sec. HVCl) . 

Alternative relaxation equation: writing the r.h.s. of 
Eq. (|D20p in the form of the divergence of a current in 
order to conserve the circulation, and using a MEPP, we 
obtain a relaxation of the form 13511 : 



dcu 
'di 



— + u • Vw = V 



D 



1 



n2{t) 



VW + l3{t)Vi) 



where ^2{t) is given by Eq. (|D29p and /3(t) by 
/ DVuj ■ VV^ dr 



m = 



n2it) J D{Vi}}fdr' 



,(D34) 



(D35) 



The boundary conditions are (^7(4)^'^ + /3(^) ^"0) • n = 
on the domain boundary. This relaxation equation satis- 
fies the same general properties as Eq. (jD28p . 



3. Relaxation equations associated with the 
maximization problem (|46p 

We shall introduce a set of relaxation equations associ- 
ated with the maximization problem (j46p . We write the 
dynamical equation as 



— +u-VuJ = X, 
at 



(D36) 



where X is an unknown quantity to be chosen so as to 
increase 5'[ZU] while conserving E and F. The time varia- 
tions of S are given by 



S = - J uJX dr. (D37) 
On the other hand, the time variations of E and F are 
E = J Xi;dr = 0, (D38) 

f = j Xdr^O. (D39) 

Following the Maximum Entropy Production Princi- 
ple, we maximize S with E ~ T = and the additional 
constraint 

?<C(r,t). (D40) 
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The variational principle can be written in the form 

SS - l3{t)6E - a{t)Sr - 



1 / X^' 



(D41) 
(D42) 



and we obtain 

X ^~D{P{t)i; + a{t)+uJ). 
Substituting Eq. ()D42p in Eq. ()D36p . we obtain 

-D{p{t)tlj + a(t)+uj), (D43) 



did „ 
— + U-WZJ 

at 



with Lo ~ ~ct{t) on the domain boundary. The Lagrange 
multipliers (3{t) and a{t) evolve so as to satisfy the con- 
straints. Substituting Eq. (|D42p in Eqs. (ID38I) and 
(jD39[) . we obtain the algebraic equations 



(D44) 



J DVuj ■ V^J dr 



(D48) 



The boundary conditions are {WuJ+/3{t)Wip)-n = on the 
domain boundary. This relaxation equation satisfies the 
same general properties as Eq. (|D43|) . If we assume that 
D is constant, the foregoing equation can be rewritten 



duj 
'di 



+ \i-Vuj = D (AZU - (}{t)uj) , 



(D49) 



2E 



E '' 



(D50) 



(V)/3(t)+«(t) = -r. (D45) 

Substituting uJ taken from Eq. (|D42p in Eq. (|D37|) and 
using the constraints (jD38[) - (|D39p . we easily obtain 

5 - J^dr, (D46) 

so that S > provided that D is positive. On the other 
hand 5' = iff X = leading to the condition of equilib- 
rium (j53p . From Lyapunov's direct method, we conclude 
that these relaxation equations tend to a maximum of en- 
tropy (or a minimum of enstrophy) at fixed energy and 
circulation. 

Alternative relaxation equation: writing the r.h.s. of 
Eq. (|D36p in the form of the divergence of a current in 
order to conserve the circulation, and using a MEPP, we 
obtain a relaxation of the form [35[: 

^ + u • VaJ = V • [Z? (VaJ + /3(t)V^)] , (D47) 



where we have used an integration by parts to obtain the 
second term of Eq. (jPSOp . 

Remark 1: Since the relaxation equations derived in 
this section solve ((TT|) . they can also be used as a numeri- 
cal algorithm to construct nonlinearly dynamically stable 
stationary solutions of the 2D Euler equations character- 
ized by a linear ui — ip relationship (see Sees. HIl and UlI)) 
independently of the statistical mechanics interpretation. 

Remark 2: Since the EHT thermodynamical equilib- 
rium with a Gaussian prior is equivalent to (|46p . the re- 
laxation equations derived in this section coincide with a 
particular case of the relaxation equations derived in [32j | 
(the ones corresponding to a Gaussian prior). 

Remark 3: Since the optimization problems ([T6| . ([36|) 
and (|46p are equivalent, the corresponding relaxation 
equations derived in Appendices ID 11 ID 21 and ID 31 have 
the same equilibrium states. However, the dynamics 
leading to these equilibrium states is different in each 
case because the constraints are different. 
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